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Abstract 

The QED radiative corrections to virtual Compton scattering (reaction 
ep — > ep^y) are calculated to first order in Oem = e^/A-K. A detailed study 
is presented for the one-loop virtual corrections and for the first order soft- 
photon emission contributions. Furthermore, a full numerical calculation is 
given for the radiative tail, corresponding with photon emission processes, 
where the photon energy is not very small compared with the lepton mo- 
menta. We compare our results with existing works on elastic electron-proton 
scattering, and show for the ep — > ep^ reaction how the observables are mod- 
ified due to these first order QED radiative corrections. We show results for 
both unpolarized and polarized observables of the virtual Compton scattering 
in the low energy region (where one is sensitive to the generalized polarizabil- 
ities of the nucleon), as well as for the deeply virtual Compton scattering. 
PACS : 13.40.Ks, 13.60.-r, 13.60.Fz 
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I. INTRODUCTION 



Virtual Compton scattering (VCS) has become in recent years a new and versatile tool 
in the study of nucleon structure and has triggered an important activity on both the the- 
oretical and experimental side (see e.g. [|l],0]). VCS, which is accessed through the (e, e'7) 
reaction, is studied now in various kinematical domains. 

At low energy, below pion production threshold, it allows to access generalized polariz- 
abilities of the proton These response functions, which constitute new nucleon structure 
observables, have been calculated in various approaches and models I^J^-^. To extract this 
nucleon structure information from VCS below pion production threshold, a considerable 
experimental effort is taking place at various electron laboratories. The first few events of 
VCS were observed in The first dedicated VCS experiment has been performed at 

MAMI and for the first time, two combinations of generalized polarizabilities have been 
determined at a four-momentum squared = 0.33 GeV^ [|TT|. An experiment at higher 



(1-2 GeV^) at JLab [0 has already been performed, which is under analysis at the time 
of writing, and a further experiment at lower is planned at MIT-Bates . 

The VCS is also studied vigorously in the Bjorken regime (where the photon virtuality 
and the photon-proton cm. energy ^/s are both large, with Q'^/s finite), which is referred 
to as deeply virtual Compton scattering (DVCS). In this region, the DVCS amplitude is 
parametrized at leading order in Q in terms of six generalized parton distributions ||I^-[16|, 
commonly denoted as skewed parton distributions (4 quark helicity conserving functions and 
2 which involve a quark helicity flip). This field has generated by now a whole theoretical 
industry, and first experiments of DVCS and related hard electroproduction reactions are 
being performed, analyzed or planned both at JLab 0, HERMES/HERA [|lD, and COM- 
PASS 

The first absolute measurement of the VCS cross section on the nucleon performed at 



MAMI [|TT|, indicates that QED radiative corrections provide an important contribution 
to the ep — >■ ep7 reaction (of the order of 20% in the kinematics considered in |11[). The 
ep — >■ ep7 reaction is particular in comparison with other electron scattering reactions be- 
cause the photon can be emitted from both the proton side (this is the VCS process which 
contains the nucleon structure information of interest) or can be emitted from one of the 
electrons (which is the parasitic Bethe-Heitler process). The radiative corrections to the 
Bethe-Heitler process are formally different compared with the case of electron scattering. 
The importance of a very good understanding of the radiative corrections is indispensable 
if one wants to extract nucleon structure information from the ep epy reaction, espe- 
cially in those kinematical situations where the Bethe-Heitler process is not negligible. The 
calculation of these QED radiative corrections to the ep epy reaction to first order in 
ctem = e^/47r ^ 1/137.036, is the subject of this paper. 

Radiative corrections were first calculated by Schwinger for potential scattering [ED 



Tsai extented the calculation of the radiative corrections to electron-proton scattering. 
The field has a long history and we refer to the standard review papers [^,^, which were 
used in the interpretation of many electron scattering experiments. 

The outline of the present paper is as follows. In section ||, we introduce the kinematics 
and notations used for the ep — >• ep'j reaction, and give the lowest order amplitudes. 

In section |T|, we give the first order QED radiative corrections to the ep epj reaction. 
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We first calculate, in section [III A| , the one-loop virtual radiative corrections originating from 
the lepton side, to the ep ejry reaction. Our strategy used to evaluate the rather com- 
plicated loop integrals, is to solve first simpler loop integrals, which contain entirely the 
ultraviolet (UV) and infrared (IR) divergences, and in which the lowest order amplitudes 
factorize. These simpler loop integrals are evaluated analytically. The finite remainder with 
respect to the original amplitude, is then expressed through Feynman parameter integrals, 
which are calculated numerically in this work. 

In section piIB| , we calculate the soft photon emission contributions from the lepton side, 
to the ep —>■ epj reaction. We discuss the similarities and differences with the bremsstrahlung 
contribution to elastic electron-nucleon scattering. These bremsstrahlung processes contain 
IR divergences which are shown to cancel exactly the IR divergences from the virtual photon 
processes. 

In section |III Q , the numerical method to evaluate the remaining finite Feynman param- 
eter integrals is presented. We discuss subsequently the cases where the integrand is regular 
or singular, the latter originating from the propagation of on-shell intermediate states in 
the one-loop corrections to the ep — > epy reaction. In particular, we discuss the different 
numerical checks performed and the accuracy of the calculation. 

In section [111 D| , we discuss the radiative corrections at the proton side and the two- 
photon exchange corrections, by referring to the recent work of Maximon and Tjon pl |. 

In section we give a full numerical calculation for elastic electron-proton scattering of 
the photon emission processes where the photon energy is not very small compared with the 
lepton momenta, and which makes up the radiative tail. We compare this full calculation 
with an approximate procedure based on the angular peaking approximation, and show to 
what extent the full calculation validates the approximate method for the case of elastic 
electron-nucleon scattering. The approximate method will be seen to be realistic enough to 
apply it next to the calculation of the radiative tail in the case of VCS. 

In section |V|, we start by briefly discussing the radiative corrections to elastic electron- 
proton scattering. We apply the radiative corrections to elastic scattering data on the proton. 
We next give our results for the ep — >■ epj reaction, and indicate how the observables are 
modified due to the first order QED radiative corrections. We discuss first the polarizability 
region for the ep —>■ epy reaction, corresponding with a low outgoing photon energy. We 
show results for both unpolarized and polarized cross sections in MAMI and JLab kine- 
matics. Subsequently, we give the effect of the first order QED radiative corrections to the 
DVCS cross section and the electron single spin asymmetry. 

Finally, we give our conclusions in section 

We present technical details needed in the calculations, in two appendices. In appendix 
we calculate the radiative corrections to elastic lepton-nucleon scattering, which serves as a 
point of comparison with the ep — >■ epy reaction. In particular, we present the details of the 
calculation of the soft photon emission contributions, and perform analytically the phase 
space integral over the soft photon in an exact way. We compare with other calculations in 
the literature. In appendix 0, we present some technical details on the integration method 
used to evaluate singular Feynman parameter integrals. 
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II. LOWEST ORDER AMPLITUDES OF THE EP EP'j REACTION 

The lowest order (in aem), contributions to the ep — > ep7 reaction are given by the 
one-photon exchange processes. We denote in this work the four-momenta of the initial 
and final electrons by k{Ee,ke) and k'{E'^,k^)] the four-momenta of the initial and final 
protons by p{E]^,'p^) and p'{E'j^,pj^); and the four-momentum of the outgoing photon by 
qWq I) <? )• Furthermore, we denote q = k — k' = p' — p + q' and = —q^ > 0. The masses 
of the electron and proton are denoted by m and respectively. The helicities of the 
initial (final) electrons are denoted by h{h'); the spins of initial (final) protons by Sp(Sp); 
and the polarization four- vector of the outgoing photon by e. The spinors of initial and final 
electrons are denoted by u{k, h) and u{k' , h'); whereas the spinors of initial and final protons 
are denoted by N{p, Sp) and N{p', Sp). Throughout this work, we follow the conventions of 
Bjorken and Drell ||25|| . 

In Figs. ^ (BHi) and (BHf), which are known as the Bethe-Heitler (BH) diagrams, a 
photon is emitted by either the incident or final electrons. The expressions for Figs. ^ (BHi) 
and (BHf) are respectively given by : 

M],^ = ze^w(fc',/z')7^i^^^^^ ^*u{k,h)j^^N{p,Sp)T,(p,p)N{p,Sp), (1) 
= ie'u{}^,h') jt* ^^^^t'^'^ Yu{k,h) iV(p,4)r,(p,p) N{p,sp), (2) 

where the electron charge is given by (— e) (i.e. e > in this work). The on-shell electro- 
magnetic vertex at the hadron side Fj, in Eqs. (^0) is given by 

F. {p',p) = F, [{p - pf) 7. + F, ((p - pf) ^^J^^^ , (3) 

where Fi and F2 are respectively the Dirac and Pauli electromagnetic (on-shell) form factors 
of the nucleon. The four-momentum squared of the virtual photon in the BH processes is 
if: = {p' ~ pYi in contrast to g^, which is the four-momentum squared for the VCS process 
7*p — »• 7p, where the final photon is emitted from the hadron side. This latter part contains 
the nucleon structure information. 

The amplitude of the VCS contribution to the e^p e~p7 reaction is given by : 

Mvcs = ~te^ u{k' ,h')^,u{k,h) \ e* H^-" . (4) 

Remark that for a positive lepton, the VCS amplitude changes sign. In Eq. (^), the gauge- 
invariant, hadronic tensor H^'^ is defined by : 

H'"' = j d'^xe-''^-^ < p'\T [f{x),f{())] \p > , (5) 

where T represents the time ordering, and j'^ the electromagnetic current operator. 

For the DVCS process in the Bjorken limit, the hadronic tensor of Eq. (^ is parametrized 
in terms of six leading twist skewed parton distributions (see e.g. [0). 
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For the VCS process at low energy, as investigated experimentally in [|Tl|-[T3[], an im- 
portant contribution to the tensor of Eq.(^ originates from the nucleon pole contributions 
shown in Figs. (BORNi) and (BORNf). The contributions of the Born diagrams to the 
hadronic tensor are given by : 

KoRN, = NiP, i) {p',p- q) ^^'^^^/""^ [P - ^^P) N(P^ ^p)' (6) 
H'soRNj = Nip, i) T'^ {p,P + q) {p + iV(p, (7) 

where the vertex F'^ is now evaluated for off mass-shell values of one of its arguments. In 
Ref. 0, the Born diagrams were evaluated by using the vertex of Eq. (^. Doing so, the Born 
diagrams are separately gauge invariant. All nucleon structure effects are then absorbed in 
a non-Born amplitude which is regular in q' and for which the Low Energy Theorem (LET) 
tells that it starts at order q' . The nucleon structure effects to the VCS tensor (Eq. (|^)) 
below pion threshold, are then parametrized at order q' in terms of six generalized (i.e. 
dependent) nucleon polarizabilities @,|]. 



III. FIRST ORDER RADIATIVE CORRECTIONS TO THE EP EP-f 

REACTION 

A. Virtual radiative corrections 

In this section, we calculate the one-loop QED virtual radiative corrections to the ep — *• 
ep7 reaction, which are represented in Fig. |^. In the present section, we consider only the 
corrections originating from the electron side as they can be calculated model-independently. 
The corrections originating from the hadronic side, for which a nucleon structure model is 
needed, will be discussed and estimated in section p.11 D| . 

The virtual radiative corrections to the BH process contain vertex corrections : Figs. |^ 
(Vli - V3i) and (Vlf - V3f); electron self-energy corrections : Figs. |^ (Si, Sf); and vacuum 
polarization corrections : Figs. ^ (Pli, Plf). We indicate in our notation of the different 
diagrams whether the photon in the ep epf reaction is emitted from the initial (i) electron 
or from the final (f) electron. 

The part of the virtual radiative corrections to the VCS process (i.e. where the pho- 
ton in the reaction ep epy is emitted from the hadronic side) which can be calculated 
model-independently, consists of the vertex diagram shown in Fig. Q (V4) and the vacuum 
polarization diagram shown in Fig. ^ (P2). The blob in those figures represents the VCS 
process. For VCS below pion threshold, the blob is given by the Born diagrams (Fig. |l| 
(BORNi) and (BORNf)) -|- non-Born diagrams, which describe the nucleon polarizability 
effects. For DVCS, the blob is given in leading order by the so-called handbag diagrams, 
where the photon hits a quark in the proton |T^,|TB|,0] . 



The calculation of the virtual radiative corrections to the VCS process is similar to 
that for electron scattering. The virtual radiative corrections to the Bethe-Heitler process 
are different, but involve the same one-loop building blocks, i.e. electron vertex, electron 
self-energy and photon self-energy. Therefore, we give in appendix ^ (sections |A 1| - [A4|) 
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the derivation and the expressions for these basic building blocks, and we apply it to elas- 
tic electron-nucleon scattering. In our calculations, we use the dimensional regularization 
method to treat both ultraviolet (UV) and infrared (IR) divergences. This amounts to eval- 
uate all loop integrals in D dimensions. The divergences then show up (when one takes 
D ^ 4) as poles of the form l/e, where e = 2 — D /2. UV divergences are regularized by 
taking D < 4 (i.e. euv = 2 — D/2 > 0), whereas IR divergences are regularized by taking 
D > A (i.e. eiR = 2 — D/2 < 0). Care has to be taken as to isolate the UV and IR divergent 
parts in the loop integrals first, as two different limits are understood when one takes D = 4 
at the end. The technical details of our calculational method can also be found in appendix 
^ We apply it here to calculate the diagrams of Fig. |^ to the ep — > ep7 reaction. 



1. Vertex correction diagrams of Figs. |^ (Vli) and (Vlf) 
The amplitude corresponding to Fig. |^ (Vli) is given by 

= (7^ (^''^) '"^ 

^ ' ^' -2k.q' ^ J (27r)^ [l^] [P - 2Lk] [P - 21. {k - q') - 2k.q'] ^ ' ^' ^ 

where a mass scale fi (renormalization scale) is introduced when passing to D 4 dimensions 
in order to keep the coupling constant dimensionless. One sees by inspection that the loop 
integral in Eq. (^, when taking D = 4, is IR finite (Z behavior), but has an UV 
divergence (/ ^ oo behavior). Our strategy to evaluate a complicated loop integral as in 
Eq. (11), is to solve first a simpler loop integral which contains entirely the UV divergence 
and which can be done analytically more easily. We observe from Eq. (H) that only the term 
in the numerator proportional to 7" /f /t* /L^a is responsible for the UV divergence. To 
evaluate it, we add a similar term by replacing one factor in the denominator and evaluate 
this term analytically. In order to obtain the equivalence with M*]^, we have to subtract the 
added term again from the expression of Eq. (H) . This leads to 



xm(A; , /i') 7 
+2 



{]k- /+m) / /■ <Pl 7° //* H 



-2k.q Y J (27r)^ \P\ \P - 2l.k] [P - 2l.k'] 



(27r)' [P] [P - 2l.k] [P - 21. {k - q) - 2k.q] 



J (27r)' [P] [P - 2l.k] [P - 2l.k'] [P - 21. {k - q') - 2k.q'] ^ ^ ' ' ' ^ ' 

It should be remarked that only the added term in Eq. (|^) (first term within curly brackets 
of Eq. d^) ) is UV divergent and has therefore to be evaluated in D dimensions using the 
dimensional regularization method. The third term within the curly brackets in Eq. (RI) is 
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the difference between the term proportional to 7" / Jt^a in Eq. (Q) and the added term 
(which has one different factor in the denominator). As can be seen by power counting, this 
term is UV finite and can therefore readily be evaluated for D = 4. The denominator in 
the UV divergent first term of Eq. @, was chosen so that it corresponds with the vertex 
correction which appears in electron scattering. Therefore, this UV divergent term can be 
calculated analytically along a similar way as was performed in appendix 0. The result is 
given by : 



4-D 



7" / ho 



(4^ 



7r 



1 



where v is defined as 



In 



- 7s + In 



47r/i2 



1 — f In 



f + 1 



)^ ^ ]k, -\- ^ ]k, -\- 



17 + 1 
V - 1 



Am? 



(10) 



(11) 



The UV divergence in Eq. ([T0|) is removed by the corresponding vertex counterterm as given 

'47r/i2' 



by Eqs. (|A^,|AT6|) 



(CT) 



« — /\/f« 

VI — ^^^BH , . x2 



(4vr) 



lE + In 
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+ 2 







- 7£; + In 



(12) 



where we have used the expression of Eq. (|l]) for the BH amplitude M^^^. Adding the 
counterterm of Eq. (0) to Eq. (P) and introducing a Feynman parametrization in the second 
and third terms of Eq. (^) in order to perform the integrals over /, yields the total, UV finite 
result : 



.5 1 



(4vr)^ 



7s + In 

^IR 



47r/i^ 



3 — f In 



v + l 
V - 1 



+ 



(4vr)^ (p - p)^ 



iV(p',Sp)r^ (p ,p) N{p,Sp) 



xu{k' , h') 7' 



(/;- /i + m) 



X 



1 



-1 



-2A;.g' 
2u 



In 



v + l 



+ 4m (£*.A;) {I - y) - jk jt* ( jk- A - 
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— 4 / dx3 x\ I dx2 X2 / dxi 



A' 
1 

with the four-vector defined by 



pi 



k - q] {I - x^) + {k- qxi) x^ x-^ 



k.q 

u{k,h), (13) 

(14) 

(15) 
(16) 



and the scalars A* and defined by 

A' = 2k.q'{l-X3) + {Pf , 

and 

B\ = 2k.q x{l -y) + m^y . 

Remark that although Eq. (|13|) is UV finite, it contains now an IR divergence through the 
vertex counterterm of Eq. (|12]) as shown in appendix A (Eq. (|A16|) ). We will demonstrate 
however in section |111 B| , that all IR divergences, arising from the one-loop corrections to 
the ep —>■ ep-f reaction, are cancelled when adding the corresponding soft photon emission 
contributions. 

The Feynman parameter integrals in Eq. ([T3| ) which orginate from the finite integrals 
in Eq. (P) remain to be evaluated. As an analytical calculation of these integrals is rather 
complicated, we will evaluate them numerically in this paper, which will be discussed in 
section Q 

In a completely similar way as for Fig. Q (Vli), the total amplitude including the coun- 
terterm corresponding to Fig. |] (Vlf) yields : 



Ml, + {CTYy, = Mi, 



(4vr) 



£lR 



- 7£; + In 



47r/i^ 



V In 



v + 1 
V - 1 



te 



(47r) {p - p) 
1 



v + l ' 



2v 



V - 1 



^/*^ + -ln 

V 



v + l 
v-l 



Jki ^ Jki ~\~ Jki ^ Jki 



2 dy dx — 7 







y (/+ 4* (/^'+ A'x) +y (X+ A'x^ 4* H 

+ 4m [e*.]^^ (l-y) - fl) H 

111 
+ 4 y dx^xl j dx2X2 j rfxi (^-^ + jpf ,t* /P^ [P^.{q-q 



((^- /i') r / 



k .q 



(17) 
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with the four- vector defined by 

P-^ = (^k' + q'^ (1 - X3) + (jt' + qxi^ X2X3 , 

and the scalars A-^ and B( defined by 

Af = -2k' .q {1 - X3) + {Pff , 

and 

B( = —2k .q x{l — y) + m?y ■ 



(18) 

(19) 
(20) 



2. Vertex correction diagrams of Figs. |^ (V2i) and (V2f) 
The amphtude corresponding to Fig. |^ (V2i) is given by 

{p -p) ^ ^ 



xu{k , h') /i 



^*u{k,h). (21) 



{27Tf [P] [l^ - 2l.k'] [P - 21. {k - q) - 2k.q] -2k.q 

One sees from Eq. (^) that again only the term in the numerator proportional to 7" 
contains an UV divergence for D = A. To evaluate the loop integral of Eq. (0), we therefore 
apply a similar trick as used before in Eq. (^. This amounts to adding and subtracting a 
term in Eq. (pl|) by replacing {P — 21. (k — q') — 2k.q') in the denominator by {P — 2l.k'), and 
which contains entirely the UV divergence. The further steps are then analogous to those 
following Eq. (^), and yield the following result for Fig. ^ (V2i) : 



(47r) 
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X u{k! , h') 
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1 



-1 + 



i,„(i±i)),,.,,i,„(i±i){,,^,V,S-4} 

2 j dy j dx^ y /l^l'' (y/i + (y/l- a') a;) + y + A') a;) Y H 
00 ^ 

+ 4m (1 -y) - {^]k- A') Y H - m V 

111 
4 y dx^xl J dx,X2 j dx, (^-^^Y + ^ rY {k - P' 

'Y r + rY > 



}(^^/^.(..) 



(22) 
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where is given as in Eq. (|T5D and where 



B2 = rn^y + x'^y (^q — + 2x k.q + 2xy k .(^q — q^ . 



(23) 



In an analogous way, the amphtude corresponding to Fig. |^ (V2f ) can be calculated, and 
yields as result : 
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1 1 



+ 4 / (ixs X3 / (ix2 X2 / dxi 



-y^ + 



where is given as in Eq. (|T9D and where 



i?2 = m'^y + x^y (^q — q^ — 2xk .q — 2xy k.(^q — q^ 



(24) 



(25) 



3. Vertex correction diagrams of Figs. |^ (V3i) and (V3f) 
The amplitude M^g corresponding to Fig. || (V3i) is given by 

mU = o N{p,s^) r. iv(p, sp) 



V3 — . , x2 



X m(/c ,h) 11 



(27r)^ [/^ + 2l.k'] [P + 2/.fc] [Z2 + 21. {k - q) - 2k.q] 



u{k,h). (26) 



Remark that the loop integral in Eq. ( PB| ) is UV finite but contains an IR divergence 
for D = 4. This is because in Fig. ^ (V3i), a soft virtual photon (/ 0) couples 
to two on-shell electron lines. To isolate the IR divergence, we first decompose the 
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numerator in Eq. (^61) by using the relations u{k\ h') (^/c' + rnj = u{k\h')2k'^ and 
{/c + m) 'yau{k, h) = 2kau{k,h). This yields : 



V3 



{P - P) 

xu{k\h')ij^-'' 



{2Tif [l^] + 2l.k'] [P + 2l.k] [P + 21. {k - q) - 2k.q] 
X |4(A;.fc') Y / + m) + A{k.k') 7" / + '^Y {/^- A' + / + rr?j ^* / /i' 

+2 ]k lY [h- /l+ /+ + Yh" {/^- /+ /* /7„} uik, h). (27) 

In Eq. (p7D, only the term in the numerator which is /-independent (the first term within the 
curly brackets) contains an IR divergence, whereas all the other terms are finite. As before, 
instead of aiming at an analytical formula for a rather complicated integral, we evaluate the 
IR divergent part of the integral in Eq. (27) by adding and subtracting a term that contains 
the divergence and that can be performed analytically rather easily. In constructing this 
term, we are looking for a denominator which contains the same dependence as the basic 
BH process in order that this BH amplitude can be factored from this IR divergent term. 
This yields the following expression, which is by construction identical to Eq. ( p7D : 

,5 




N{p,s [p,p) N{p,Sp) 



X 



A{k.k')Y {A- / +m) 
^ [P] [P + 2l.k'] [P + 2l.k] [-2k.q] 

1 

(27r)^ [P] [P + 2/.A;'] [P + 2l.k] [P + 21. {k - q) - 2k.q 

P -21. {k-q) 



A{k.k')Y [jk- +m) ^ 



-2k.q 



+ A{k.k')Y / 



+ 7" /7''(^-^'+/ + m) /7a]} w(A;,/i). (28) 

Remark that the added term (first term of Eq. (^)) contains the IR divergence whereas the 
other terms of Eq. (^8])) do not have any divergences so that the corresponding integrals 
may be performed directly in four dimensions as indicated. For the first term of Eq. ( P5| ) we 
furthermore see that the /-independent part of the energy denominator is the same as the 
one occuring in the corresponding Bethe-Heitler diagram (Fig. |l]a). The /-dependent part 
of the energy denominator for this term is the same as the one for the vertex correction to 
elastic electron scattering, Eq. (|A9|) . The corresponding integral may therefore be evaluated 
analytically in a similar way as was done in appendix A. This yields for the IR divergent 
term in Eq. (|28| ) : 
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■ (29) 



The evaluation of the finite four- dimensional integral in Eq. (^) can be performed at the 
expense of the introduction of three Feynman parameter integrals due to the four energy 
denominators : 

1 



/21 \12 



[P + 2l.k'] [P + 2l.k] [P + 21. {k - q) - 2k.q] 
11 1 



6 dyy^ dx2 X2 / dx 



1 



with the four-vector -P^^^j defined by 

^xiX2 = [q - xi q'^ X2 + k' , 
and the scalar C* defined by 

^ 2k.qx,X2 + y {P^,,S . 
The final result for the amplitude Myg is then given by 



4' 



(30) 



(31) 



(32) 
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(33) 
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The Feynman parameter integrals in Eq. (|33D will be performed numerically as explained in 
Section pTCl 



In an analogous way, the result for the amplitude My^ corresponding to Fig. ^ (V3f) can 
be calculated, and yields as result : 



Ml, = Mij, 



(4vr) 



— 7£; + In ( 
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2v 
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V - 1 
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2k'. q 



2 



1 



(cfy 



Cf (Cf) 

+ 1/ (-aT'^ (/+ /*+ 4m5-) /^ij} 

with the four- vector P^^^^ defined by 

PLx2 = - {q - xiq'^ X2 + k , 
and the scalar defined by 

= -2k'.qxiX2 + y (p4J' . 



(34) 



(35) 



(36) 



Remark that in the vertex correction diagrams where the photon couples to the final 
electron (diagrams of Fig. |^ denoted by /), the invariant mass of the virtual (e" + 7*) 
state in the loop is given by + 2k .q > m? . This means that an on-shell propagation 
is possible for the (e~ + 7*) state. This translates mathematically into the presence of 
integrable singularities in the corresponding Feynman parameter integrals of Eqs. (0, ^4D 
and (P^), and yields an imaginary part for the corresponding amplitude. In contrast, in 
the vertex correction diagrams where the photon couples to the initial electron (diagrams of 
Fig. ^ denoted by i), the invariant mass of the virtual (e~ + 7*) system in the loop is given 
by which means that the corresponding integrals contain no singularities. The numerical 



treatment of those singular Feynman parameter integrals will be discussed in section [III C 
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4- Electron self-energy diagrams of Figs. ^ (Si) and (Sf) 



We next evaluate the electron self-energy diagrams of Figs. (Si) and (Sf). We only have 
to consider those diagrams where a photon is emitted and re-absorbed by an intermediate 
electron line. The diagrams with a loop on the initial or final electron lines are already 
absorbed in the wavefunction and electron mass renormalization, and therefore do not yield 
an additional correction. This can also be seen from the expression Eq. ( |A30 ) for the 
renormalized lepton self-energy, which vanishes on-shell. 

The amplitude corresponding to Fig. |^ (Si) is then given by : 

Msi = ie^ , ,^ ,2 ^{p. Sp) r. (p,p) N{p, sp) 



(37) 



(p - pY 

X n(fc,/i07" ^ , , ' S(A;-gO— — , , ' ^*u{k,h), 



-2k.q 



-2k.q 



where the renormalized self-energy is denoted by S and is given by Eq. ( [A3 01 ). Remark 
that the UV divergence in the loop integral of Fig. ^ (Si) has been removed through the 
renormalization of the electron field and electron mass. The UV finite renormalized self- 
energy E contains however an IR divergence from the counterterms. Inserting the expression 
for t (Eq. ( |A30|) ) into Eq.(|3), yields : 

^2 r 1 /4^^2- 
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The amplitude corresponding to Fig. |^ (Sf) is given by : 
1 
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which can be worked out analogously as before and yields : 
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Note that in Fig. ^ (Sf), the four-momentum squared of the (e~ + 7*) state in the loop 
is given by {k' + g')^ = + 2k'. q' > m?. Therefore, the self-energy and the amplitude 
for Fig. (Sf) is complex, as was also noted for the vertex diagrams of Fig. ^ where the 
photon is emitted from the final electron (denoted by /). Eq. (^0]) yields indeed a complex 
amplitude because In (— 2/c'.g'/m^) = \Yi{2k' .q' /rn?) + ztt, for k' .q' > 0. 



5. Vertex correction diagram of Fig. |^ (V4) 

The vertex correction to the VCS process is given by Fig. | (V4), and its calculation is 
the same as the one for elastic electron scattering. This yields for the renormalized vertex 
correction : 



Mv4 = -ie^u{k',h') 



(F(Q2) _ f(Q2 = 0)) 7. - GiQ^) la,^ 



q'^ 
2m 



u{k,h)-elH^r (41) 



In Eq. (^TJ), F{Q^) — F^Q"^ = 0) is given by Eq. ( |A17|) and reduces in the ultrarelativistic 
limit (Q^ >> m"^) to Eq. ( |A19|) . The magnetic correction G{Q'^) is given by Eq. ( All ), and 
vanishes in the ultrarelativistic limit. 



6. Vacuum polarization diagrams of Figs. |^ (PH, Plf) and (P2) 

The vacuum polarization corrections of Figs. |^ (Pli, Plf) and (P2) involve the renor- 
malized photon self-energy n((5^), which has been calculated in appendix Therefore, we 
get for the vacuum polarization correction to the BH process (Figs. |^ (Pli, Plf)) : 

ML = MIh i , and = M^^ 1 , (42) 

^1 ^^i_n(-t)' ^^i_n(-t)' ^ ^ 

with t = {p' — pY- Similarly, we get for the vacuum polarization correction to the VCS 
process (Fig. ^ (P2)) : 

Mp2 = Mvcs -^TT^iT • (43) 

1 - n(Q^) 

In the ultrarelativistic limit (Q^ >> m^), II(Q^) is obtained from Eq. (|A40|) 



B. Soft-photon emission contributions and cancellation of IR divergences 

After removing the UV divergences from the virtual photon corrections to the ep — > eprf 
reaction in the last section, the resulting expressions still contain IR divergences. Both the 
corrections to the BH process of Figs. | (Vli, Vlf, V2i, V2f, V3i, V3f, Si and Sf ) and the 
vertex correction of Fig. |^ (V4) to the VCS process contain IR divergences. It is known 
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for QED since a long time p6| , |27[] , that these IR divergences are cancelled at the cross 
section level by soft photon emission contributions. These soft photons are emitted from the 
charged particle lines and can have energies up to some maximal value AEg which is related 
to the finite resolution of the detector. In appendix (section |A 5| ), we calculate the soft 
bremsstrahlung contribution to electron scattering by performing the phase space integral 
over the soft photon in an exact way, and give the finite correction (after cancellation of all 
IR divergences) to the elastic electron scattering cross section. In this section, we generalize 
the result of appendix ^ to the case of the ep — >■ epj reaction. The diagrams for the 
ep — > ep~f reaction with one additional soft photon are shown in Fig. ^, where the hard 
photon of the ep — > epy process is indicated by its four-momentum q'. In this section, we 
will show that the soft photon emission contributions of Fig. ^ contain IR divergences which 
exactly cancel the IR divergences appearing in the virtual photon correction diagrams of 
Fig. 1^. The process where the energy AEg of the additionally emitted photon is not very 
small compared with the lepton momenta in the process, makes up the radiative tail to the 
ep — i> ep7 reaction. Its calculation will be discussed in section [TV. 



1. Factorization of amplitude for soft-photon emission processes 



Here, we evaluate the diagrams of Fig. ^ in the soft photon limit, i.e. when the second 
emitted photon has an energy much smaller than the initial and final lepton energies and 
also smaller than the hard photon (denoted by q') in order to distinguish both photons. We 
will see that only the diagrams where a soft photon couples to an on-shell lepton contain IR 
divergences and lead to a finite logarithmic correction in AEg. 

The amplitude corresponding with Fig. ^ (bli) is given by : 



Mf 



bli 



ie^u{k', h') 7 



. (A-^-/+^) 
-2k.q' -2l.{k - q') 



,^ {jk- j[ + m) 
-2k.l 



-e /*(/)) u{k,h) 



X 



{p - pY 



N{p ,Sj,)T^ [p ,p] N{p 



'1 "pj ) 



(45) 



where / is the four-momentum of the soft photon. In the soft photon limit (/ — > 0), 

(2A;°- /7°) u{k,h) ^ 2k''u{k,h), 

(46) 



Eq. (^51) simplifies by using (/;— / + m) 7" m(A;, /i) 
which yields for Eq. (|i5| ) in the soft photon limit : 



M, 



bli 



M' 



BH 



-e)el{l) 



where M]j^ is the Bethe-Heitler amplitude of Eq. (^ - corresponding with photon emission 
from the initial lepton. Similarly, we can derive the amplitude for Figs. |^ (b2i, blf and b2f) 
which yields in the soft photon limit : 



Mbif + Mb2f 
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BH 



BH 



-e)4(0 
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k'" 
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(47) 
(48) 
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where M^fj is the Bethe-Heitler amphtude of Eq. - corresponding with photon emission 
from the final lepton. 

Figs. ^ (b3i) and (b3f) contain the contributions where the soft photon couples to an 
off-shell lepton line. The amplitude corresponding with Fig. ^ (b3i) is given by : 

l\/r ■ 4-/;' ii\ V Kh" ti~ l + m) /l' + to) I 

^ 2 N{p', 4) Nip, sp) , (49) 



(P - P) 

In the soft photon limit, Eq. (^9]) can be simplified by using : 



-2k.q' -2l.{k-q')^ -2k.q' ~ -2k.q' ^ -2k.q' 

{jk~ /i' + m) {k~q'Y Y 



(50) 



-2k.q' ~k.q' -2k.q' 
Consequently, the amplitude of Eq. (^) is given by : 

+ ie^u{k' ,h') '^ ^-2kq' ^'^^ ^^^'^^ ^^"i ^^'P ' ^'v^^'' N{p,Sp). (51) 
Similarly, the amplitude corresponding with Fig. |^ (b3f) is given by : 

+ I e'uik\ h') ^*^'^'lfj^^ u{k, h) ^-^^ N{p, .;) (p,p) Nip,sp). (52) 

In complete analogy to Eqs. ( ^^8]) , we can also calculate the soft photon emission 
contributions to the VCS process. They are shown in Figs. ^ (b4) and (b5), and their 
calculation in the soft photon limit yields : 



Mb4 + Mb5 = Mvcs{-e)eUl) 



¥l~kl 



(53) 



where Mycs is the VCS amplitude of Eq. (^). 

We see from Eqs. ( ^ , ^ and ^), that for the diagrams of Fig. ^ where the soft photon 
couples to an on-shell lepton, the original amplitude factorizes : in Eqs. (|^,|^) the BH 
amplitude factorizes, and in Eq. (^) the VCS amplitude factorizes. The resulting amplitudes 
are proportional to 1//, which leads to a logarithmic divergence when integrating over the 
phase space of the soft photon. In contrast, the amplitudes of Eqs. (|5T|j5^) where the photon 
couples to an off-shell lepton line are finite when Z — 0, and the corresponding phase space 
integral becomes vanishingly small in the limit 1^0. 
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2. Radiative correction due to soft-photon emission processes 



In the soft-photon hmit, we therefore need only to keep the bremsstrahlung corrections 
of Eqs. p7| , |48| and ^3]), where the BH and VCS amphtudes factorize. To first order in ojem 
(relative to the BH + VCS cross section) the bremsstrahlung correction therefore amounts 
to calculate the phase space integral of the form : 



da 
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(54) 



where 1 = |/| denotes the soft photon energy, and where the total BH amplitude is given by 
Mbh = M^BH + ^BH- '^^^ calculation of the bremsstrahlung integral of Eq. (|^) goes along 
similar lines as the corresponding integral for elastic scattering, for which the technical 
details can be found in appendix ^ (section |A5D . We will point out in this section the 
differences which arise for the ep —>■ epj reaction. 

There are two practical ways to measure the ep —>■ epj reaction, by measuring two 
particles in the final state. One can either measure the outgoing electron in coincidence 
with the recoiling nucleon : this is the ideal technique when measuring the ep —>■ ep^ 
reaction at low outgoing photon energy as is done in pAj-pTSf . The alternative is to measure 
the outgoing electron in coincidence with the photon : this is the technique when doing a 
very inelastic experiment, such as deeply virtual Compton scattering, where the photon is 
produced with a large energy. We discuss here first the case where one detects the outgoing 
electron and photon, and indicate at the end the changes which apply when measuring the 
outgoing electron and recoiling nucleon. 

If one measures the ep —>■ epj reaction by detecting the outgoing electron and photon, 
one eliminates in Eq. ( [A41|) the integral over with the momentum conserving (5-function, 
which gives : 



da 



d'k' 



d^q 



dH 



1 



(27r)' 2E'^{2'kY 2\q\ (27r)' 21 2E'^ 



N 



X {271)5 [Ee + En - E'^-\q 



X \Mbh + Mvcs? (-e') 



1 1 

\(1 \ - 




PN - q - ly + Ml 




kfj. 








k'.l 


' kl 




¥1 " kl_ 





(55) 



Due to the energy conserving 5- function in Eq. ( |A42| ) , the upper limit in the integration over 
the soft photon phase space depends on the angle. Therefore, this integration volume has 
a complicated ellipsoidal shape in the lab system. In order for the soft-photon phase space 
integration volume to be spherical, one has to perform the calculation in the cm. system Si 
of the (recoiling nucleon + soft-photon), generalizing the procedure of appendix ^ for elastic 
scattering to the ep — > epry reaction. The system Si is defined by : p^ + l = Pn + Q — Q = 0. 
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In the system Si, the energy conserving delta function in Eq. (|55|) is independent of the soft- 
photon angles, and the maximal soft photon energy is isotropic. The integral over the soft- 
photon momentum (up to some maximum value AEg) can then be performed independently 
from the integration over the soft photon emission angles. If AEg is sufficiently small, one 
can furthermore neglect the soft photon energy with respect to the other energies in the 5- 
function, and perform the integral over the photon momentum |g | in Eq. (^) to obtain the 
correction to the fivefold differential ep — > epy cross section. We indicate in the following 
only how the squared matrix element for the ep —>■ epj reaction is modified due to soft 
photon emission. This correction due to soft bremsstrahlung is given by : 



r d^i 








1 {27i f 21 






¥1 ~ kl_ 



(56) 



^ gives immediately the correction factor 
, the soft-photon phase space integral is 



The correction factor multiplying \Mbh + Mycs 
to the fivefold ep epj cross section. In Eq. 
understood to be performed in the system Si, where the integration volume is spherical. Its 
calculation was already performed in appendix 0. One sees that the integral in Eq. (^) has 
a logarithmic IR divergence, corresponding with the emission of photons with zero energy. 
To evaluate it, one has to regularize it, which is done in this work by using dimensional 
regularization. This amounts to evaluate the integral (in the system Si) in D — 1 dimensions 
{D —>■ A corresponds to the physical limit). This calculation is performed in appendix ^ and 
yields (similar to Eq. ( [A63D ) as result : 



vcs\ 
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^IR 
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(57) 



In Eq. (pTf), 5r is the finite part of the real radiative correction corresponding with soft 
photon emission, and is given as in appendix ^ (Eq. ( [A66| )) by : 



tr 




71 

y 



+ Sp cos^ 




(58) 



In Eq. (|58|), we next have to express the kinematical variables {E^., E'^,cos6e) in the 
system Si (denoted by tilded quantities) in terms of the lab quantities, which we denote by 
untilded quantities {Ee, E'^, cos9e)- To make the transformation between the system Si and 
the lab system, we first introduce the missing four-momentum p^i = Pjy + l. The system Si 
is defined by Pmi = 0, and the missing mass M^i of the system {p' + I) is defined by : 



Mi 



ml 



ip' + i? = {v + q- q'f 



(59) 



We can then easily express the electron energies and angle in the system Si in terms of lab 
quantities : 
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Ee = ^^ = ^k.{p + q-q) = -—\E,-— — , (60) 



sin^ ej2 = S& sin^ eJ2 . (62) 
The maximal soft-photon energy AEg in the system Si, is given by : 

If one measures the ep —>■ epy reaction by detecting the outgoing electron and recoiling 
proton, the derivation goes along similar lines as above. One starts now by eliminating in 
Eq. ( [A41| ) the integral over q . Then one goes into the cm. system ^2 of the (VCS photon 
q' + soft photon), where the energy conserving 5- function is independent of the soft-photon 
angles, and where the maximal soft photon energy is isotropic. This system S2 is defined 
by: q + I = p^ + q — Pj^f = 0. The calculation of the soft-photon emission integral is 
then completely similar as above, and leads to the finite correction of Eq . (^8|) , where the 
the kinematical variables {E^, E'^, cosOe) are now understood in the system (S2. To make the 
transformation between the system ^2 and the lab system, we first introduce the missing 
four-momentum pm2 = Q + The system ^2 is defined by Pm2 = 0, and the missing mass 
Mm2 of the system (g' -|- I) is defined by : 

Ml2 = iQ' + if = ip + q-pT- (64) 

We can then easily express the electron energies and angle in the system ^2 in terms of lab 
quantities : 
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, ^ ^ ^A;'.(p + g-p') = Mil(e'^ + -^- ^\ , (66) 



sin2 ~ej2 = sin^ ej2 . (67) 



The maximal soft-photon energy AE^ in the system 1S2, is given by : 



AE, = ^, (68) 



3. Cancellation of IR divergences 

We can now demonstrate for the ep — > ep7 reaction, that the IR divergences from the 
soft photon emission corrections exactly cancel against the IR divergences from the virtual 
radiative corrections, calculated in section [III A| . Concentrating here only on the IR divergent 
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parts of the virtual radiative corrections, we found in section [III A] that the amphtudes of 
Eqs. (p!3| , |T7| , ^ , p^ , |33| , |3^ , |38| , ^ , and ^l]) contain IR divergences. Those IR divergent parts are 
given by : 
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(69) 
(70) 
(71) 
(72) 
(73) 



Adding them all up gives the following correction to the squared amplitude for the virtual 
radiative corrections : 



= |Mbb + M 



vcs 





- 1 



7b + In 



+ ..., (74) 



where the ellipses denote the finite first order virtual radiative correction to the ep — >■ ep7 
reaction as was calculated and can be found in section [111 A| . Adding the virtual (Eq.([7^)) 
and real (Eq.(|57D) radiative corrections to the ep eprf reaction, one verifies that the IR 
divergences in the sum exactly cancel, showing QED at work! Note that this cancellation 
is different as compared to the case of elastic electron scattering. Indeed, for the virtual 
photon correction diagrams to the Bethe Heitler process, there are 3 types of vertex dia- 

r2f)), and the corresponding 



grams (Eqs. ( p^JTOj and ^Tj)) and the self energy diagram (Eq. 
counterterms, which have an IR divergence. On the other hand, for the virtual radiative 
corrections to elastic electron scattering, there is only one vertex diagram which is IR diver- 
gent. 



C. Integration method for the virtual photon corrections 

At this stage of the calculation of the first order QED radiative corrections to the 
ep — i> eprf reaction, the treatment of all UV and IR divergences, resulting from the ra- 
diative corrections at the electron side, has been performed. The UV divergences have been 
removed by the renormalization procedure whereas the IR divergences were shown to cancel 
at the cross section level when adding the soft photon emission processes. Now, the evalu- 
ation of the remaining Feynman parameter integrals in the finite terms such as in Eq. ( p!3D 
still has to be done. 

Among the one-loop virtual radiative corrections to the ep — > epy reaction shown in 
Fig. 0, six give rise to simple analytical formulas. For the six vertex diagrams, denoted by 
Vli,V2i,V3i,Vlf,V2f and V3f, the trick consisting of adding and subtracting the divergent 
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term for each of them (as explained in Section |III A| ) gives rise to Feynman parameter inte- 
grals that are rather complicated to be done analytically. Therefore, we will evaluate them 
in this work by a numerical procedure. Although these Feynman parameter integrals are by 
construction finite, appropriate numerical methods are needed to perform them. Two main 
difficulties are encountered in these numerical integrations. Firstly, the variations of the 
integrated functions are always extremely sharp near the integration limits. In fact, a typi- 
cal behaviour is a rather fiat dependence in the middle of the domain and two pronounced 
rises when approaching or 1 for the Feynman parameters with a width of the order m/ E^.. 
The contribution of these two peaks has to be evaluated carefully in order to obtain a good 
precision for the final result. Secondly, we know that the virtual radiative corrections to the 
ep — i> ep7 reaction allow the propagation of on-shell states (see section [III A|) . This is math- 
ematically expressed by the presence of integrable singularities in the Feynman parameter 
integrals which require an analytical continuation into the complex plane and gives rise to 
an imaginary part for the amplitude. 

To evaluate the Feynman parameter integrals, our strategy is to perform the first in- 
tegration analytically. The last integrations will then be performed numerically using the 
Gauss-Legendre integration routine. The analytical calculation of the first integration pro- 
vides a shorter calculational time and a higher precision. The main advantage however is 
that in the case of a singularity, the pole is avoided by deforming the integration contour into 
the complex plane, using analytical continuation. In this way, one removes the difficulties 
for the remaining integrations along the real axis. 

To classify the Feynman parameter integrals that occur in the six vertex diagrams under 
study, we start by factorizing all the Dirac 7 matrices and decomposing the components of 
the four-vectors. All resulting integrals then reduce to the generic form : 

1 

dxidx2dx3 1^^' ^ ^75") 

Q{Xi,X2,X3) 



where P and Q are polynomials in three (real) Feynman parameters Xi,X2,X3. Let's choose 
Xi to be the more internal variable. Then the first integration is either of the form : 
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xT dx^ ^^g^ 
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xT dxi ^^^^ 







{axl + + 7)" 



where a, (3 and 7 are polynomials in X2 and 0:3 with coefficients that are functions of 
kinematical variables. In Eqs. (^,^), m varies from to 4 and n is equal to 1 or 2, to 
accomodate all cases appearing in section [III A| . These successive decompositions increase 



the number of terms to calculate but they have the advantage to provide two simple classes of 
integrals without any vector or matrix dependence. The possibility of poles in the integrands 



of Eqs. ([76[ , [77[) naturally splits the problem into two parts, whether the integrand is regular 
or singular. 
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1. Regular integrand 



When the denominator doesn't have any singularities, some recurrence relations exist for 



these integrals and can be found in Ref. |^8[. Unfortunately for small values of a as compared 
to (3 or to 7, it has been seen that these relations are numerically unstable. This has thus 
led us to use several methods of integration with each a different domain of validity. For 
small ratios r {r = a/ (5 for Eq. (|76D or r = 0/7 for Eq. (|77D ) as compared to 1, we perform 
a Taylor expansion of the integral and tune the order of each development to complete a 
fixed criterium of convergence (for example we require that the ratio between the last and 
the first terms is of the order of the numerical precision in double precision). For r > 1 the 
recurrence relations are used as they are stable in this range. In the intermediate zone 
(0.2 < r < 1), we use the Gauss-Legendre numerical integration method. 



2. Singular integrand 

In the case of the propagation of on-shell intermediate states, the polynomials of the 
denominators in Eqs. (|7^,^) acquire one (or two) roots in the domain of integration. 
Some simple physical considerations have shown that among the six diagrams numerically 
evaluated, the three processes where the photon in the ep — > ep7 reaction is emitted from 
the initial electron line are free of poles (section |III A| ). In contrast, the three vertex graphs 



where the photon is emitted from the final electron line were seen to contain singularities. 
The corresponding integrals are then defined by an analytical continuation into the complex 
plane and take the form : 

1 1 

dX\ f Xi dX\ frrr,\ 







[axi + /? ± ie)"- J {axi + fixi + 7 ± 



The prescription for on-shell propagation is of course already taken into account in the 
propagators and determines the sign in front of the ie (which can also be obtained by 
applying the simple trick m m — ie/2). Complications can occur from the possibility of 
two distinct roots in the interval [0, 1] for the second order polynomial. An important remark 
then concerns the variable of integration. In Eq. (|75D , the choice of xi as the more internal 
dimension was purely arbitrary. In fact all the decompositions in the three parameters 
have been derived and it has been shown that it was always possible to find an expansion 
providing at most one singularity. 

In appendix P, we give the analytical results for the integrals of Eq. (^). We checked 
these results with a numerical method, where one pole along the interval [0,1] is avoided by 
analytically continuing the integrand into the complex plane. In this way, the integral along 
[0,1] is replaced by an integration along a semi-circle (with origin at 0.5 + 0i and radius 1/2) 
in the opposite complex half-plane with respect to the pole. A comparison between the two 
methods shows a perfect agreement. Only in the special cases where one pole comes close 
to an edge of the domain of integration [0,1] (typically within a distance m? / to or 1), 
one has to increase the number of integration points of the numerical method to obtain the 
same precision. 
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3. Numerical checks and accuracy 

Thanks to the analytical calculation of the first integration in the Feynman parameter 
integrals under study, singularities on the real axis have been removed and the two remaining 
integrations can then be performed numerically using the Gauss-Legendre method. In the 
implementation of this algorithm the major difficulty consisted in finding the suitable binnig 
of the integration domain and to determine the number of points per bins. A detailed study 
of the integrated functions has been performed to estimate the width and amplitude of the 
sharp variations close to the ends of the domain. In this paragraph we discuss various checks 
of the precision of our results as well as their numerical stability. 

A strong cross check of the reliability of our calculations is the exact agreement between 
two programs developed in paralell Both of them use the same numerical method 

but they have been coded independently using in most cases different decomposition of 
the terms and different order in the integration variables, which checks the symmetry in 
the permutation of Xi,X2 and 0:3 variables. Comparison at each intermediate stage of the 
calculation also excludes any missprints in the writing of the quite extensive expressions. 

Besides this agreement between two independent programs, the next requirement is the 
numerical convergence of the calculations. Figs. § and ^ show results obtained for typical 
MAMI and JLab kinematics respectively. Beyond a certain density of integration bins and 
points per bin, the numerical instabilities are brought down to few 10~^ of the lowest order 
cross section. This accuracy is far below all the other theoretical uncertainties related to the 
performed approximations or experimental knowledge of the form factors (of the order of 
1%). Nevertheless this kind of very good convergence is useful since numerical instabilities 
can be amplified in the coherent sum of all the diagrams or when computing higher energy 
kinematics. In the case of the deeply virtual Compton scattering, we have checked that one 
has to double the number of integration points, to get the same numerical precision. 

Some features of the electromagnetic interaction itself can also be used to check further 
the validity of our results. Let's consider the total amplitude of the sum of all the virtual 
radiative correction diagrams (Fig. Denoting the Lorentz index associated with the real 
photon vertex by /x, this amplitude can be written as the scalar product T'^e* where e* stands 
for the polarization vector of the real photon with four-vector q and where represents 
the electromagnetic current. The gauge invariance of electromagnetism implies T^q^ = 
and provides us with a powerful test of our calculations. Since our numerical accuracy is 
finite, we cannot get exactly zero. Therefore, we rather define a quantity compared to which 
the scalar product T'^q^ has to be small. A natural quantity is the product of the norms of 
the two Lorentz vectors. The gauge invariance criteria thus becomes a test of the smallness 
of the following dimensionless ratio : 





2 




(g'0)2 



This ratio is shown in Fig. ^as a function of the angle between q and q . The gauge invariance 
is verified by the fact that the smallest ratio (solid curve) stays in the range [10~^, 10~^] and 
is obtained when the complete set of diagrams with analytical+numerical terms is included 
in T^. 
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As a last consistency check, we investigated the mass dependence of the virtual radiative 
corrections. The relative effect in the BH + Born cross section is illustrated in Fig. ^ for 
different values of the mass of the lepton. For this test we kept track of the mass dependence 
in all the kinematical variables. We observe that when increasing the lepton mass (at fixed 
lepton kinematics), the effect of the radiative corrections rapidly decreases, which reflects 
the suppression of photon emission by a heavy particle. 



D. Radiative corrections at the proton side and two-photon exchange corrections 
In section [111 A| - pi C[ we calculated the radiative corrections to the ep —>■ epj reaction. 



corresponding with the diagrams of Figs. |^ and ^ They are the virtual radiative correc- 
tions at the lepton side, the vacuum polarization contributions and the soft-photon emission 
from the lepton. These can be calculated model-independently as has been shown above. 
Although these corrections are the dominant ones (when >> m^, leading to large loga- 
rithms), we want to estimate in this section how large are the virtual radiative corrections 
at the proton side, the two-photon exchange corrections (direct and crossed box diagrams) 
and the soft-photon emission from the proton. Generally, the radiative corrections from the 
proton side are typically suppressed compared with those from the electron, due to the much 
larger mass of the proton. However, to calculate the first order radiative corrections to the 
ep — >■ epj reaction which originate from the proton side, one needs a model for the VCS 
process. We do not aim in this paper, to calculate these corrections within a given model. 



However, to provide some quantitative estimate, we will follow the results of , where the 
corrections at the proton side were studied for elastic scattering. 

The Z-dependent corrections originate from the interference between soft-photon emis- 
sion from the electron and from the proton, and from the two-photon exchange contributions 
(direct and crossed box diagrams). Both processes contain IR divergences, which cancel in 
their sum at the cross section level. The interference between the soft-photon emission from 
the electron and from the proton can be calculated along the same lines as in appendix |A 5| 
for the electron (neglecting form factor effects in the soft-photon limit). For the two-photon 
exchange contributions, the calculation is dominated by those regions in the integration 
where one of the two exchanged photons is soft. Therefore, one can evaluate the rest of this 
amplitude by taking the momentum of either of the two exchanged photons to be zero. In 
this approximation, the original amplitude factorizes and one can follow the derivation of 
P^ , where this same calculation has been performed for elastic scattering. Therefore, the 
Z-dependent radiative corrections can be estimated in the soft-photon limit by the same 
correction factor of Eq. ( |A75D as for elastic scattering. 



The Z^-dependent corrections originate from the soft-bremsstrahlung from the proton 
and from the proton vertex corrections. In ||24], these corrections have also been calculated 



for elastic scattering. For the soft-photon emission, one can again factorize the orginal am- 
plitude, so that the same correction factor is obtained for the ep — epy reaction as for 



elastic scattering. The proton vertex correction has been split in |^ into two parts. The 
first part contains entirely the IR divergence, which cancels with the IR divergence from 
soft-photon emission from the proton, and in which the original amplitude factorizes. The 
second term in the proton vertex correction depends on the nucleon structure (form factor 
dependence for elastic scattering) and will be different when going from elastic scattering to 
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the ep — >■ ep7 reaction. For elastic scattering, this structure dependent term was however 
found [24] to be quite small, except when going to very large (much larger than Mf^). 
When staying in the few GeV^ region, this correction was calculated in 
1%. Therefore, we approximate the dependent correction to the ep 
the structure-independent term of Eq. ( |A77|) , as calculated in 
following the structure dependent term. 



to be well below 
ep7 reaction by 
and will neglect in the 



IV. RADIATIVE TAIL FOR ELASTIC SCATTERING AND VCS 

Besides the knowledge of the virtual radiative corrections and the soft-photon emission 
contributions to the ep — > ep7 reaction, which were studied in section |Tl|, the accurate de- 
termination of the ep —>■ ep7 cross section from measured spectra also implies the knowledge 
of the radiative tail. The radiative tail consists of the photon emission processes where a 
semi-hard photon (with energy not very small compared with e.g. the lepton energies) is 
radiated from the electron (or proton). 

The radiative tail to elastic or inelastic lepton-nucleon scattering has been the subject 



of numerous studies in the literature 31,p2lE4,R3i. The elastic radiative tail also makes a 



sizeable contribution to the cross sections for deep-inelastic lepton-nucleon scattering (see 
e.g. H). 

One should notice that the distinction between the soft-photon emission and the radia- 
tive tail is not a fundamental one, the latter being just the extension of photon emission 
processes to higher energies. Although the formulas given in this paper for the real radiative 
corrections can in principle be extended and applied to higher energies (e.g. Eqs. ( |A71[ - 
A73| ) for elastic scattering), in some cases the characteristics of the experimental detection 



apparatus can be such that the cut in E'f'' — E'^ (elastic case) or in the missing mass Mj 
(M^]^ or for the VCS case) cannot be cleanly defined, because the apparatus can have 
a changing acceptance as a function of E'^^ — E'^ or Mj, introducing a bias in the radiative 
tail. Therefore, it is useful to consider the radiative tail separately and to generate it in 
a Monte Carlo simulation. In doing such a simulation, it can be very helpful to have a 
"recipe", because it is a way to fold radiative effects with acceptance functions and other 
effects (e.g. multiple scattering, energy loss by collision, external radiative effects). In the 
literature such "recipes" were quite often presented. Many of them are based on one or 
another version of the peaking approximation, introduced originally by Schijf . In the 
peaking approximation, the photon is radiated along either the initial or final electron di- 
rections, i.e. the direction of the electron is not changed while radiating, only its energy is 
changed. 

Below, we start by giving such a recipe, based on the formulas presented in this paper. 
What one essentially needs for a Monte Carlo simulation is an electron energy loss distri- 
bution due to real internal radiative effects. For each event one can then sample in such 
a distribution, both for the incoming and the outgoing electron. We next give a compari- 
son between such a method based on the peaking approximation, with an exact numerical 
calculation of the radiative tail. We show to what extent the full calculation validates the 
approximate method for the case of elastic electron-nucleon scattering, and show that this 
method is realistic enough to apply it next to the calculation of the radiative tail in the case 
of the VCS. 
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A. Energy loss distribution for real internal radiative effects 



The details of the calculation of the real radiative corrections can be found in appendix 0. 
It is discussed there how the real internal radiative corrections give rise to a correction factor 
e^^ to the cross section. The part of 6r giving rise to the radiative tail (when differentiating 
6 a with respect to the electron energy loss) is the first term of Eq. ( |A66| ) , which contains 
the maximal energy of the emitted photon AEg, which is defined as in Eq. ( A49| ) . The 
correction factor e*^^ can be written as the product of a number of factors, of which the first 
one is given by : 



2\ a 



EeK 



(80) 



where a is given by (see Eq. ( |A66| )) : 



TT 



.ni«; 



(81) 



and where the tilded quantities in Eq. (^) are expressed in the cm. system of (soft pho- 
ton + recoiling proton) as explained in appendix |A 5| . Because in a simulation it is more 
straightforward to apply radiative effects in the lab, we express Eq. ( pO]) in lab quantities, 
by using Eq. ( |A49|) , which yields : 
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where AE'^ = E'^^ — E'^. Introducing furthermore the quantity AEe = rfAE'^, we can write 
Eq. (PD as : 
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(83) 



The energy changes AE^ {AE'^) can be interpreted as the energy losses of the incoming 
(outgoing) electron due to radiation before (after) the scattering process respectively. We 

can then interpret the factor ^-^^ j as the fraction of incoming electrons which have lost 

an energy between and AEe, after being subject to real internal radiation in an equiv- 

f AE' \ " 

alent radiator with thickness a. The factor ( -^t^ ) has a similar interpretation, but then 

on the outgoing electron side Q. Given this interpretation, if one uses a AE distribution 
Iint{E , AE , a) , which satisfies : 



Note that when applying Eq. (83) to the radiative tail, i.e. when considering the emission of 
a photon whose energy is not very small compared with the electron energies, we calculate E'J''' 
in the formula for AE'^, using the elastic scattered energy corresponding with an initial electron 
which has radiated and whose energy is given by Eg — AE^. In the soft-photon limit this difference 
disappears. 
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I,r,t{E,AE,a)d{AE)= l^j , (84) 

then it is clear that by samphng such a distribution in a Monte Carlo simulation, the 
correction factor is correctly obtained. The distribution I^nt, which has this property is 
given by : 

/,„i(E,AE,a) = ^(^-^j , (85) 

and is normalized to 1 : 

\ I,nt{E,AE,a)d{AE) = l. (86) 
Jo 



B. Evaluation of the radiative tail and comparison with an exact numerical 
calculation for elastic electron-proton scattering 

Given the above distribution, a method for introducing a radiative tail due to internal 
radiation in a Monte Carlo simulation for elastic electron scattering suggests itself : 

i) For the incoming electron, sample an energy loss AE^ using the distribution (^) with 
E = Ef. the incoming electron energy. 

ii) Apply elastic electron scattering using the reduced electron energy E^, — AE^,, and if 
the cross section behavior is implemented in the simulation, use the elastic scattering cross 
section at the reduced electron energy. After the elastic scattering process, the outgoing 
electron has an energy E'f''. 

iii) For the outgoing electron, sample an energy loss AE'^ using the distribution ( p5D with 
E = E'f. The final electron energy is now E'^^ — AE'^. 

To calculate the equivalent radiator thickness a of Eq. (|8TD, one needs the value of 
Q^, which one can in principle only calculate after the complete process has taken place. 
However, one can show that the above procedure reproduces the correction factor (^) with 
a very good accuracy already by calculating the value of with elastic electron scattering 
kinematics. 

It is intuitively clear that the above procedure, in the case where a constant cross section 
is used, will reproduce the correction factor of Eq. (^). In case the actual elastic scattering 
cross section behavior is implemented, the cross section "walk" with the incoming electron 
energy is taken into account. Remark that the above procedure implies an electron energy 
loss both at the incoming and the outgoing electron sides. 

The discussed method implies, however, the assumption of a strict alignment of the 
bremsstrahlung photons in the direction of the radiating leptons, which is known as the 
(angular) peaking approximation. The strength on the other hand is found by integrating the 
correct angular shape in the soft photon limit, as done in appendix [A 5| . To test the validity 



of this approximate procedure, we performed a fully numerical calculation of the radiative 
tail for elastic electron-proton scattering. It consists of integrating over the photon phase 
space in the diagrams where a photon is emitted from an electron (cfr. BH diagrams of Fig. |l| 
(a) and (b)), as well as the diagrams where a photon is emitted from the nucleon (cfr. Born 
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diagrams of Fig. |l] (c) and (d)). In doing so, we nowhere neglect the photon momentum /, in 
contrast to the calculation of appendix A_5 in the soft-photon limit, where this momentum 
is neglected with respect to the lepton momenta. For fixed electron kinematics, the angular 
phase space of the soft photon is covered by a grid with about 225000 points, chosen with 
increased density in the peak regions in order to keep the point-to-point change of the cross 
section smaller than 10 %. Attention has to be paid right in the middle of the peaks where 
the cross section drops very rapidly to (practically) zero within the characteristic angle 



m/Ee, as shown in Fig. |[ More details on this numerical integration can be found in []36 
The result of this integration is the absolute cross section of the radiative tail, differential 
in the outgoing electron's momentum and angles. It is shown by the points in Fig. ^ for 
Ee = 855.0 MeV and 6e = 52.18°. The energy of the outgoing electron is then determined 
by E'^ = E'^^ — AE'^. The points are compared with the analytical result in the soft-photon 
limit, obtained by differentiating the expression of Eq. ( |A66|) for 6r - for photon emission 
from the electron - with respect to AE'^. This gives a strict AE^T^ behaviour, yielding the 
cross section cTq = aBomd/ AE'^ where the proportionality factor a is given as in Eq. (^TJ). 
The soft-photon formula gives thus a straight line when both the cross section and AE'^ 
are presented on a logarithmic scale. The deviation can be seen in the lower plot of Fig. |^. 
From the keV-region up to about 1 MeV for AE'^, the deviation is less then 10~^ which 
can be taken as an upper limit for the error of the numerical integration procedure. This 
agreement demonstrates that the soft-photon approximation holds to very good precision in 
this region. For higher values of AE'^, a raise of the photon emission cross section is observed 
as is expected due to the change of kinematics leading to a lower momentum transfer to the 
proton, and to a resulting "walk" of the cross section. We also show on the lower plot of 
Fig. ^ the result when both radiation from the electron and proton are considered. For better 
presentation, both results are normalized to the cross section 0"^ for soft-photon emission 
from the electron, as defined above. 

In Fig. 10, we compare for two kinematics the exact numerical calculation of the radiative 
tail with the approximate method of the Monte Carlo simulation as discussed above. The 
simulation has been investigated by running it with and without the cross section behaviour 
(dipole form factors assumed), and the ratio between the two versions is presented by the 
lines, the outer lines representing the statistcal accuracy. One notices that the increase 
of the radiative tail is reproduced, but somewhat overestimated compared with the exact 
calculation. 



C. Application to virtual Compton scattering 

The above procedure can also be applied to VCS, as long as the angular peaking ap- 
proximation is used, i.e. the electron does not change its direction while losing energy by 
internal real radiation. Indeed, Eq. ( p8D is completely similar to the elastic case, when ex- 
pressing it in the cm. system of either (soft photon -|- outgoing nucleon) or (soft-photon 
-|- outgoing photon) depending on how the ep — > ep7 reaction is measured, as explained 
in section [III B| . After exponentiation, one can apply a factorization completely similar as 
in Eq. (PB|). Because under the assumption of the angular peaking approximation AEg/E 
is constant under a Lorentz transformation, we obtain the property that the shape of the 
distribution ( |S3| ) is system independent, only its endpoint value E changes. As a result. 
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one can apply the distribution of Eq. (|85|) in the lab for VCS, but then using lab values 
for Ef, and E'^. For VCS, one certainly can have a changing acceptance of the detection 
apparatus as a function of missing mass (making a "clean" cut in missing mass on the data 
impossible), so that generating a radiative tail in a Monte Carlo simulation with the above 
described method is probably the best way to implement the radiative tail correction to the 
data. Such a simulation was implemented for the VCS experiments already performed at 
MAMI [O] and at JLab |T2|, and will be fully described in a forthcoming paper 



V. RESULTS AND DISCUSSION 

A. Elastic electron-proton scattering 

Before showing results for VCS, we briefly discuss first the effect of the radiative correc- 
tions to elastic electron-proton scattering, in order to have a point of reference. The radia- 
tive corrections to elastic electron-proton scattering are presented in detail in appendix 0. 
In Table |, we show for different elastic kinematics (MAMI, JLab) the numerical values 
of the vertex correction {S^ertex of Eq. (|A68| )), the vacuum polarization correction (5„ac of 
Eq. ( |A70|) ), and the real radiative correction at the electron side {Sr of Eq. ( |A66| )). We also 
show the Z and dependent corrections, 6i (Eq. (|A75[ )) and 6^^ (Eq. ( |A771 )) respectively, 
as derived in the recent work of [^. We omit here the small part in the dependent 
correction which depends on the particular model for the nucleon structure (in the elastic 
case, the form factors), as can be found in [Q. In Table |, we indicate the total radiative 



correction 6tot as the sum of all the different contributions as in Eq. ([A74| ) . From Table |, we 



see that by far the largest contribution to the radiative correction comes from the large log- 
arithm and double logarithm in Q"^ /rn? in the electron vertex correction. When evaluating 
the real radiative corrections for E'^^ — E'^ = 0.01 i?e, the total effect of the radiative correc- 
tion is an upwards correction of the data (for negative 5tot) of the order 20 - 25 %. In the 
last column of Table | (denoted by EXP), we also indicate the result when exponentiating 
all corrections except the vacuum polarization contribution, which - as modification of the 
photon propagator - is resummed as in Eq. ( A73| ) . One sees that this can lead to differences 



of the order of 2 %. 

In an elastic scattering experiment, one measures an scattered electron spectrum and 
has to evaluate the real radiative corrections as a function of the cut {E'^^ — E'^ which one 
performs in the spectrum. Dividing the measured cross section by the correction factor 
(1 -|- 5tot) and plotting the result as function of {E'^^ — E'^), should then lead to a "plateau" 
behavior, which demonstrates the consistency of the procedure (within a certain range of 
the value {E'^'' — E'^) where one knows the radiative tail to sufficient accuracy). 

The determination of the elastic cross section for the kinematics Ep = 705.11 MeV, 6^ 



= 42.6° is shown in Fig. 0. The upper plot shows the dE-spectrum of elastic data taken 
(during the beam time of the VCS experiment) at MAMI. It is compared with the simulated 
spectrum (dashed line). On the lower plot, the ratio of the experimental spectrum integrated 
up to the value AE'^, to the simulation integrated also up to AE'^ is shown as function of the 
cutoff energy AE'^. This gives the elastic cross section, which is seen to be stable below the 
1%-level over a long interval up to the cut by the acceptance of the spectrometer. The slow 
descent for higher AE'^ indicates that the simulation overestimates slightly the radiative tail. 
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B. VCS below pion production threshold 



We next turn to the ep — >■ ep7 reaction below pion threshold. It was discussed in section 
|I|, that the lowest order (in aem) amplitude of the ep ep^ process at low outgoing photon 
energies q' = |g | is given by the BH + Born processes. The deviation from the BH + Born 
amplitudes grows with q', and can be parametrized (at low q') in terms of six generalized 
polarizabilities (GP's) of the nucleon, which are function of . 

A first VCS experiment has been performed at MAMI [|1T|. It consisted of measuring 
the ep — i> ep7 reaction at five values of q' below pion threshold, ranging from q' = 33 MeV/c 
to q' = 111.5 MeV/c. At the lowest value q' = 33 MeV/c, where the polarizability effect is 
negligeably small, the measurement serves as a check of the Low Energy Theorem (LET). 
The measured deviation as function of q' can then be interpreted as the effect of the GP's. 
It is clear that both to test the LET as well as to extract the GP's from the measured 
deviation with respect to the BH + Born result (which is expected to be of the order 10 - 
20 % at the highest q' value), it is a prerequisite to know very accurately how the result is 
modified due to radiative corrections. 

In Fig. 0, we first show the differential cross section for MAMI kinematics at a low 
value q' = 33 MeV/c, as function of the cm. angle of the emitted real photon with respect 
to the direction of the virtual photon. One sees from Fig. that the virtual radiative 
corrections reduce the BH + Born result in these kinematics by about 16 % (or when 
applied to data, increase the uncorrected data by 16 %). The real radiative corrections have 
to be estimated as function of the cut which one performs in the missing mass spectrum. 
The VCS experiments below pion threshold measure the ep — *• eprf reaction by detecting the 
outgoing electron and proton, and reconstruct the missing mass as defined in Eq. (|6^. 
In Fig. |13|, the real radiative corrections are shown for a value of AE^ = 10 MeV, where 
the soft-photon energy AEg is determined from the cut in the missing mass according to 
Eq. (|68|). For the small value q' = 33 MeV/c, the real radiative correction depends only 
very little on the angle (through the last terms on the rhs of Eqs. ( |65| , |66D ) . For AEg 
= 10 MeV, the real radiative correction is given hy Sr ^ - 0.025, which corresponds 
with increasing the uncorrected data by about 2.5 %. For AEg = 20 MeV, ~ + 0.02 
(reducing the uncorrected data by about 2 %), and for AEg = 30 MeV, 6^ ^ + 0.045 
(reducing the uncorrected data by about 4.5 %). To determine the ep — *• epj cross section 
from the measured missing mass spectra, one has to perform a consistency check by plotting 
the experimentally measured (uncorrected) cross section divided by the radiative correction 
factor as function of the cut in the missing mass spectrum. In this way, one has to find a 
"plateau" behavior, as was demonstrated before for elastic data. This consistency check was 
also performed on the VCS data measured at MAMI |TT|, and will be shown in a forthcoming 
publication |]37[| . 

In Fig. |T^, we show the the differential cross section for MAMI kinematics at the highest 
measured outgoing photon energy : q' = 111.5 MeV/c. The virtual radiative corrections are 
mainly q' independent (for these rather small values) and lead thus also here to a reduction 
of the BH + Born result by about 16 %. The real radiative corrections are again shown for 
AEg = 10 MeV, and exhibit a slight angular dependence. These corrections were applied 



to the data from the unpolarized MAMI experiment of |11|. From the deviation of the 



radiatively corrected data and the BH + Born result, two combinations of GP's have been 
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extracted at - 0.33 GeV^ in |TT] 



An experiment below pion production threshold to measure the GP's at higher has 
also been measured at JLab [|12| and is under analysis at the time of writing. In Fig. |1^, we 
show how the BH + Born cross section is modified due to the virtual radiative corrections. 
It is seen that for the JLab kinematics of Fig. |T5|, the BH + Born result is reduced at the 
backward angles by about 20 % due to the virtual radiative corrections. 

The unpolarized VCS cross section below pion threshold provides three independent 
structure functions (when varying the value of e in the experiment), which allows to extract 
three of the six (lowest order) generalized nucleon polarizabilities. To extract the three 
remaining nucleon polarizabilities, one has to resort to double polarization observables as 
discussed in . In particular, double polarization observables with polarized electron beam 
and with a polarized target (along either of the three axes), or alternatively by measuring 
the recoil nucleon polarization, provide three new observables to extract the three additional 
nucleon response functions [^,0- In Fig. W, we show the double polarization asymmetries 
for MAMI kinematics, by measuring the recoil polarization components along the 2;-direction 
(virtual photon direction) or along the x-direction (perpendicular to the virtual photon but 
parallel to the scattering plane). One aims to extract the polarizability effect in these 
observables from the deviation of the measured asymmetry and the BH + Born result (see 



e.g. ||38| for an estimate of this effect within a model calculation). Therefore, it is important 
to know how much the BH + Born result is affected by the radiative corrections before 
extracting the polarizability effect. It is seen in Fig. ^ that the effect of the radiative 
corrections on the double polarization asymmetries nearly drops out in the ratio (much less 
than 1 % change of the asymmetries). At the low values of the outgoing photon energy q' 
(e.g. q' ~ 33 MeV/c) where the polarizability effect is very small, these asymmetries are 
also hardly affected by radiative corrections. Therefore, these asymmetries can also provide 
an independent check of the LET. An experiment to measure the VCS double polarization 
observables by measuring the recoil nucleon polarization is planned at MAMI in the near 
future |3l. 



C. Deeply virtual Compton scattering 

Besides the low energy region, the VCS process is also studied in the Bjorken regime, 
where and u = p.q/M]^ are large, with a;^ = Q"^/ {2M f^v) fixed. In this kinematical region, 
the process is refered to as deeply virtual Compton scattering (DVCS). In the Bjorken regime, 
the DVCS amplitude factorizes into a perturbatively calculable hard scattering amplitude, 
and into a non-perturbative part at the proton side, expressed in terms of so-called skewed 
parton distributions (SPD's) which generalize the ordinary parton distributions. These 
SPD's are new nucleon structure observables which one aims to extract by measuring e.g. 
the exclusive ep epy reaction in the Bjorken regime. Similarly as was seen before in 
the threshold region, the ep eprf reaction can have an important contribution from the 
BH process, besides the DVCS process of actual interest. However, the BH and DVCS 
contributions behave differently as function of the lepton beam energy, as studied in Refs. 
ro|,|^PT|]. In particular, at the lower beam energies, such as e.g. available at JLab, the BH 



process dominates in the forward direction over the DVCS process. In this region, the DVCS 
process becomes only measurable due to its interference with the BH process. In order to 
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extract the DVCS process (and the nucleon structure information) from its interference with 
the BH, it is therefore important to have good knowledge of how the radiative corrections 
modify the BH amphtude. 

In Fig. 0, we show the ep ejrj cross section in kinematics accessible at JLab, where 
such an experiment is planned [|T^. The DVCS cross section is calculated by using the 
ansatz for the SPD's of It is seen from Fig. [l^, that the BH indeed dominates over the 
DVCS cross section in these kinematics, and that the DVCS cross section gets enhanced due 
to its interference with the BH. One furthermore sees that the virtual radiative corrections 
reduce the BH + DVCS cross section by about 23% in these kinematics. This is mainly due 
to the reduction of the BH process when including virtual radiative corrections. The real 
radiative corrections are shown in Fig. |l^ for a value /S.Es = 0.1 GeV, which corresponds 
with a cut in the recoiling hadronic missing mass spectrum (defined in Eq. (p9D) of - 
Ml ~ 0.21 GeV2. 

In Ref. 0, it was suggested that an exploratory study of the DVCS process might be 
possible by studying the ep — ep7 reaction with a polarized electron beam. The electron 
single spin asymmetry (SSA) does not vanish out of plane and is only due to the interference 
of the BH amplitude and the imaginary part of the DVCS amplitude (i.e. the BH amplitude 
does not lead to a SSA, because it is purely real). Therefore, one expects this SSA to be 
less sensitive to radiative corrections on the BH amplitude. However, as the BH amplitude 
enters the SSA linearly in the numerator, but quadratically in the denominator (as in the 
unpolarized cross section), one might wonder what is the residual effect of the radiative 
corrections on this observable. In Fig. |l^, we show the SSA for DVCS at JLab. One sees 
that the SSA gets only slightly reduced due to the radiative corrections. The reduction of 
the SSA amounts to maximum 5% of its value around 5°, where the asymmetry reaches its 
maximal value. Therefore, the SSA shows to be a rather "clean" observable for extracting 
the DVCS amplitude in a region where the BH process dominates. Its measurement is also 
envisaged at JLab in the near future [O. 



VI. CONCLUSIONS 

We studied in this work the first order QED radiative corrections to the ep — > eprf 
reaction. The one-loop virtual radiative corrections have been evaluated by a combined 
analytical-numerical method. Several tests were shown to cross-check the numerical method 
used. Furthermore, it was shown how all IR divergences cancel when adding the soft-photon 
emission processes. A fully numerical method was presented for the photon emission pro- 
cesses where the photon energy is not very small compared with the electron energies, which 
makes up the radiative tail. Besides, we have also presented an approximate calculation of 
the radiative tail, which was shown to be realistic enough for use in a Monte Carlo simula- 
tion. 

We compared our results first to elastic electron-proton scattering. Subsequently, the 
results for the radiative corrections to the ep —>■ epj reaction were shown both below pion 
threshold and in the deeply virtual Compton scattering regime. 

Below pion threshold, our calculations were apphed to the first dedicated VCS experiment 
at MAMI, and show that the effect of the radiative corrections results in an enhancement of 
the uncorrected data by about 20 % (or an equivalent reduction of the theory). VCS double 
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polarization asymmetries where shown to be insensitive to radiative corrections. 

For the DVCS, we calculated radiative corrections for JLab kinematics and found the 
virtual radiative corrections to lead to an enhancement of the data by about 23 %. The 
single spin asymmetry was shown to be only slightly reduced by radiative corrections. 

Although we focussed here on the kinematical regimes of ongoing or planned experiments, 
the present work can also serve as a tool in the analysis of future VCS experiments. 
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APPENDIX A: RADIATIVE CORRECTIONS TO ELASTIC LEPTON-NUCLEON 
SCATTERING USING THE DIMENSIONAL REGULARIZATION METHOD 
FOR BOTH UV AND IR DIVERGENCES. 



In this appendix, we provide the reader with some details of the derivation of the radiative 
corrections to elastic lepton scattering at one-loop level. In our derivation, we use the 
dimensional regularization procedure to regularize both ultraviolet and infrared divergences. 

After a short introduction of the renormalization method, we calculate subsequently the 
vertex diagram at the lepton side (Fig. |19|(a)), the lepton self-energy diagram (Fig. |19|(b)), 
the vacuum polarization diagram (Fig. |l^(c)), and give an analytical result, without ap- 
proximations, for the soft photon emission at the lepton side (Fig. |1^ (d) and (e)). We 
compare our results with other derivations found in the literature. At the end we collect 
the results to correct the elastic lepton-nucleon scattering cross sections and discuss the role 
of the radiative corrections at the proton side and the two-photon exchange corrections by 
referring to the recent work of Ref. P^. In this appendix, we use the same notations as 
explained in section ||. 



1. Renormalization method 

In calculating QED radiative corrections in this work, we are using the BPHZ renormal- 
ization method (as explained e.g. in Ref. [0), which consists of replacing in the unrenor- 
malized Lagrangian all bare quantities by renormalized ones. For QED, the bare Lagrangian 
is given by (we are using the conventions of Bjorken and Drell in this work) 



CB = ^Bitrd^-mB)<ilB - \Fb^.Fb'"' - CB^BY^BABf., (Al) 

where the bare field tensor Fj^'^ is given by 

Fj^" = 9^ - d" A% . (A2) 

The renormalization of the theory amounts in redefining the bare quantities in terms of 
renormalized (i.e. physical observable) ones : 

= Zl'^ ^ , Ab^ = Zl'^ A^" , 

ruB = Zm m , cb = Zg e . (A3) 



In Eq. (^), the renormalized finite quantities are ^,A^,m and e. A theory in which all 
divergences can be absorbed into renormalization constants such as Z2, Z^, Zm and Zg in 
Eq. ([A3|), is called multiplicatively renormalizable. This procedure leads to a decomposition 
of the QED Lagrangian of Eq. ([Al|) into 



Ct3 = Cn + CcT, (A4) 
where C-ji represents the renormalized Lagrangian in terms of the physical (finite) quantities 

Cn = ^{t7^d^-m)<If - ^F^.F'^" - e^7^^^M> (A5) 
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and where Lcr is called the counterterm Lagrangian 

£cr = (^2 - 1) ^^t'^S^^ - {Z^Z^-\)^m^ - {Z^ - \)-F^,F^^ - {Z^-\)e^^^^A^ . 

(A6) 

In Eq. ( |A6| ) , the vertex renormalization constant Zx is defined as Z^ = ZgZ2Zl^^. For a 
renormalizable theory such as QED, all divergences obtained by calculating loop diagrams 
with the renormalized Lagrangian C-ji are cancelled by the corresponding contributions in 
the counterterm Lagrangian Car- It will be shown below how the QED renormalization 
constants are calculated to order O(e^) by calculating the vertex diagram, the lepton self- 
energy diagram and the photon polarization diagram at the one-loop level. 

As QED is a gauge invariant theory, we will simplify all calculations in this work by 
using the Feynman gauge. 



2. Vertex diagram 



The on-shell photon-lepton-lepton vertex is represented by 



= uik\h') 



leA" k,k 



u{k, h) , 



(A7) 



and the on-shell vertex of Eq. can be parametrized as 



uik! ,h') k^" (k! ,k^ u(k,h) = u(k',h') \(1 + F(Q^)) Y - G(Q^)ta^"'—] u(k,h) 
V / ^ 2mJ 



(A8) 



where q = k — k' . 

To order O(e^) , the vertex A'^ (corresponding with Fig. |T9|(a)) is given by 



k ,k 



{2n) 



D 



[P + 2l.k'] [/2 + 2l.k] 



+ O ie') , (A9) 



where a mass scale (renormalization scale) has to be introduced when passing to D ^ A 
dimensions in order to keep the coupling constant dimensionless. It is immediately seen by 
power counting that in four dimensions (D = 4), the one- loop integral in Eq. ( |A9| ) contains 
an ultraviolet (/ oo) logarithmic divergence and an infrared (/ 0) logarithmic diver- 
gence. To subtract the divergent parts (by the corresponding counterterms) of expressions 
such as Eq. ( |A9|) , one has to regularize them first. 

We follow in this work the dimensional regularization procedure to regularize both ultra- 
violet and infrared divergences. The dimensional regularization method amounts in calcu- 
lating loop diagrams in D dimensions. Physical observables are obtained by letting — ^ 4 
at the end. To obtain an integral which is ultraviolet convergent, one has to take D < 4, or 
ejjv = 2 — Z}/2>0in expressions such as Eq. ( [A9[ ). To obtain an integral which is infrared 
convergent, one has to take D > 4, or eiR = 2 — D/2 < 0. The two different limits show that 
care has to be taken with the limit D A, which means that the parts in Eq. ([A9|) that are 
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infrared divergent and the parts that are ultraviolet divergent have to be separated and in 
the corresponding terms, two different limits have to be taken when one approaches D = 4. 
Although the dimensional regularization scheme has been applied originally to ultraviolet 
divergent expressions as it respects the symmetries of the theory (in particular the gauge 
symmetry for a gauge theory), it has also been applied in a few works to regularize infrared 
divergences p3| , p |. 

When working out the integral in Eq. ( |A9D , one obtains after some algebra the following 
expressions for -F(Q^) and G{Q'^) to order O(e^) : 



(4vr) 
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where v is given by 



4m^ 



Q 



2 ' 



(AlO) 



(All) 



(A12) 



with Q"^ = - q"^ > 0. In Eq. ( |A1CI| ), je represents the Euler constant, and the Spence (or 
dilogarithmic) function is defined by 



Sp{x) 



dt 



ln(l - t) 
t 



(A13) 



From Eq. ( [ATT|) , the one-loop radiative correction to the electron magnetic moment 
follows as 



(AM) 



which is the result first obtained by Schwinger |^ 



To remove the UV divergence from the vertex correction Eq. ( |A1CI| ), one has to determine 
the vertex renormalization constant Zi of Eq. ( [A6[ ). Zi is determined by requiring that the 
total vertex 



A/^ = + (Zi - 1)7^ 
defines the physical electron charge at = 0, i.e. 



(A15) 



= 1 - F(Q2 = 0) 
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It is seen that the vertex renormahzation constant Z\ contains besides the UV divergence 
also an IR divergence. The renormahzed vertex of Eq. (|A15 ), is determined by the vertex 
correction function -F(Q^) — F{Q'^ = 0) which is given to first order in (where aem = 
e^/47r) by the expression 



F(g2)_F(g2 = o) 
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(A17) 



The expression for the vertex correction function F{Q'^) — F{Q'^ = 0), which was calcu- 
lated here using the dimensional regularization method for both the UV and IR divergences, 
agrees with the ones derived in many textbooks (see e.g. Eq. (47.52) of Ref. where a 
full derivation is given). The correspondence with the calculations which use a finite photon 
mass (A) as IR regulator is found to be 



1 , /47r^2 

eiR \ 



< > In^ . 



(A18) 



In the ultrarelativistic limit (Q^ >> m?), the vertex correction F{Q'^) — F{Q^ = 0) can 
be found from Eq. (|A17|) to be given by 



F(g2) - F(g2 = 0) 
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(A19) 



It is seen from Eq. ( |A19| ), that the finite part of the vertex correction at high is dominated 
by a quadratic logarithmic term. 



3. Lepton self-energy diagram 

The free lepton propagator (for a lepton with four-momentum k) 

S^k) ^ + "^ 



(A20) 



^2 _ ^2 _|_ 

is modified through the lepton self-energy S (k), to the full lepton propagator 

S{k) = S"{k) + S°{k) S(A;) S{k) . (A21) 
To first order O(e^), the lepton self-energy (corresponding with Fig.|T^(b)) is given by 



iT.{k) 



2 4- 

-e iT 



(27r)^ [l^] [{k + iy-w?] • 



(A22) 
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By power counting, it is seen that the integral of Eq. ([A22|) contains a hnear UV divergence 
but is IR finite in the hmit D ^ 4. The integral of Eq. ( A22 ) can be worked out and yields 
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(A23) 

where fc^ = k'^/rn^. 

To remove the UV divergence from the self-energy Eq. ([A23D , one has to determine the 
renormalization constants Z2 and from Eq. ( |A6| ). This counterterm contribution leads 
to the renormalized self-energy 



S (fc) = S {k) - {Z2-I) jk + {Z^Zrn - 1) m . 



(A24) 



Inserting Eq. ( [A24| ) into Eq. ( |A2 1| ) and developing S (/c) as a Taylor series expansion around 
jk = m yields for inverse of the total lepton propagator 



= {^-m) 



1 - 
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+ 0((/;-m)2) 



(A25) 



Requiring that the total propagator S has a pole at = m with residue 1, determines the 
renormalization constants Z^ and Z^ as 



1 + 



{l-Zm)Z2m = S(/; = m). 
Using the first order expression of Eq. (|A23|) for the lepton self-energy, yields 



(A26) 
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Remark that although the unrenormalized lepton self-energy of Eq. ([A22| ) is IR finite, 
the lepton field renormalization constant Z2 contains an infrared divergence for the derivative 
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of S that appears in its definition (see Eq. ([A271 )). Furthermore, a comparison of the first 
order expressions for the lepton field renormahzation constant Z2 (Eq. (|A28|) ) with the 
vertex renormahzation constant (Eq. ( |A16| )) shows that they are the same (It is known 
as a Ward identity and can be shown to hold to all orders as a consequence of the gauge 
invariance of QED). 

Finally, using the expressions of Eqs. (|A28|) , (|A29|) , the renormalized lepton self-energy to 
first order in a^m is given by 
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. (A30) 



It is seen from Eq. ( |A3CI|) that for an on-shell lepton {/c = m), the renormalized lepton self- 
energy E is exactly zero. Consequently, this correction has only to be applied for internal 
lepton lines. 



4. Vacuum polarization diagram 



Starting form the free propagator of a photon with four-momentum q (as stated before, 
we give all expressions in the Feynman gauge) 



q2 



the full photon propagator can be written as 

D^'^iq) = Dr(g) + D^^iq) U,,{q) D'/{q) 



(A31) 



(A32) 



where Uf^xlq) represents the vacuum polarization correction. To order 0(6^), the vacuum 
polarization (corresponding with Fig.[l9|(c)) due to I" loops (with lepton / = e,fi,T) is 
given by 



The gauge invariance of QED leads to the relation q'^q'^ Il^^x^q) = (Ward-Takahashi iden- 
tity). Consequently, the vacuum polarization correction can be written as 



n.A(g) = {-9nxq^ + q.qx) H (g^) 



(A34) 



where the function 11 (g^) is IR convergent and contains only a logarithmic UV divergence 
as can be seen from Eq. ( |A33| ). 

Using Eq. ( |A34|) , the self-consistent relation for the full photon propagator (Eq. ( |A32| )) 
yields 



D^^iq) 
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g2 (l-n(g2)) 
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+ term m q^ q , 



(A35) 



where we don't have to specify the term in q'^ q'^, as the photon propagator will be contracted 
with conserved currents on both sides so that this term in q'^ q" will not contribute to physical 
observables. Evaluating the one- loop integral of Eq. ( [A33| ) , one obtains 
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where v is given by Eq. ( |A12| ). 

The UV divergent term in Eq. ( |A36|) is removed by adding the counterterm in of 
Eq. ([A^). This leads to the renormalized photon propagator 



D^'^iq) 



+ term in q^ q^ 



g2 (^l-n(g2) 

where the renormalized photon polarization 11 is given by 

n(g2) = n(g2) _ - 1) . 



(A37) 



(A38) 



Requiring that the renormalized photon propagator (Eq. ( [A35| ) ) has a pole at g2 = g with 
residue 1, determines the renormalization constant : 



l + n(g2 = 0) . 



(A39) 



Consequently, the renormalized finite photon polarization is found from Eqs. ( |A36| ) and ( |A38| ) 
to be given by 
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which agrees with the result derived in Ref. [it . 



5. Soft photon emission contributions 



The calculation of the one-loop vertex correction of Eq. ( [A9| ) was seen to be both UV 
and IR divergent. The ultraviolet divergence was removed by renormalizing the fields and 
parameters of the theory. The remaining infrared divergences are cancelled at the cross 
section level by the soft bremsstrahlung contributions p6| , p7| . 

In this bremsstralung process (see Figs.|T^ (d) and (e)), an electron is accompanied by 
the emission of a soft photon of maximal energy Ai?s (which is related to the detector 
resolution and is therefore much smaller than the electron energy which radiates this soft 
photon). To first order in (relative to the Born cross section) the bremsstrahlung cross 
section amounts to calculate a phase space integral of the form : 
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where 1 = |/| denotes the soft photon energy, and where Mborn denotes the Born amphtude 
for elastic lepton-nucleon scattering. In Eq. ( |A41|) , terms in the soft photon momentum 
were neglected compared with the electron momenta k and k , except in the denominators 
of the lepton propagators where they matter. 

If one performs an experiment where the outgoing electron is detected, and where the 
recoiling proton remains undetected (i.e. if one measures a single arm electron spectrum), 
one eliminates in Eq. ( |A41|) the integral over with the momentum conserving 5-function, 
which gives : 
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Due the energy conserving 5-function in Eq. ( [A42|) , the integration volume for the soft 
photon has a complicated ellipsoidal shape in the lab system. In order for the soft-photon 
phase space integration volume to be spherical, one has to perform the calculation in the 
cm. system S of the (recoiling nucleon + soft-photon), as discussed in |2T[]. The system S is 
thus defined by : p^ + l = Q + Pn = 0. In the system 5, the energy conserving delta function 
is independent of the soft-photon angles, and the maximal soft photon energy is isotropic. 
The integral over the soft-photon momentum (up to some maximum value AE^) can then 
be performed independently from the integration over the soft photon emission angles. If 
AEg is sufficiently small, one can furthermore neglect the soft photon energy with respect to 
the other energies in the 5-function, and perform the integral over the electron momentum 
|/Cg I in Eq. ( |A42|) . The integration over the outgoing electron momentum eliminates the S- 
function, which implies the elastic scattering constraint. This yields then for the differential 
cross section with respect to the outgoing electron angles, the following correction due to 
soft bremsstrahlung : 



da 



da 
d^t 



BORN 



dH 



{2tiY 21 



k'.l 



kfj, 
kl 



k'f" 
¥1 



(A43) 



REAL SOFT-t 

where the soft-photon phase space integral is performed in the system S, in which the 
integration volume is spherical. We will denote in the following the external kinematics 
in the system S by tilded quantities {Ee, E'^, En, E'^^) to distinguish them from the lab 
quantities, which we denote by untilded quantities {Ef>, E'^, E^ = MN,E'jy). To make the 
transformation between the system <S and the lab system, we first introduce the missing 
four- momentum p^ = p'jy + 1. The system S is defined by Pm = 0, and the soft photon limit 
implies ~ M^- We can then easily express in the system S, the energies for the external 
particles in the elastic scattering process, in terms of lab quantities : 
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where the elastic scattering condition (Q^ = 2M7v(i?e — E'^)) has been used in the last step 
in Eqs. ( A44 , [X45| ). The angle 6e in the frame S is obtained from k.k' = E^E'^ll — cos^e) = 



EeE'^{l — cosOe), which shows (using Eqs. ( |A44|jA45D ) that in the soft-photon limit, this 



angle is the same as in the lab system, i.e. cos6'e = cos^e- 

The integral of Eq. ( |A43| ) extends up to a maximal soft-photon energy AEg in the system 
S, which is expressed in terms of the lab quantities Ee and E'^, by using : 

ip' + lf-Mlj = {p + k- k'f ~Mlj = 2p-{k- k') + {k- k'f , (A47) 

which leads (for soft-photon energies, i.e. keeping only terms of first order in AEs) to 

2MnAE, ^ 2MN{Ee - K) - AEeE'^ sin^ , 

= 2Mn{E, - E'^) - 2Mjv(^e - E'f)E'jE'f . (A48) 



All quantities on the rhs of Eq. ( |A48|) are in the lab, and the elastic scattering condition has 
been used in the last line {E'J''' denotes the elastic scattered electron lab energy, to distinguish 
it from E'^). From Eq. (|A48| ) , one determines then AEg in terms of lab quantities from the 
scattered electron spectrum through 

AEg = V {E'f - K) , (A49) 



where the recoil factor 77 is given by 77 = E^/E^ 
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Deviations from the soft-photon emission formula Eq. ( |A43| ) will show up when AEg is 



not very small compared with the lepton momenta in the process. The emission of such a 
semi-hard photon is what is usually referred to as the radiative tail. Although the distinction 
is somewhat arbitrary, one can always split the integral for photon emission into two parts, 
one by integrating up to a small value AEg, where the soft-photon approximation in writing 
down Eq. (|A43| ) holds, and a second integral, starting from this small (but non-zero) value 



of AEg up to the energy where one performs the cut in the spectrum. This second integral 
is finite and can be performed numerically. Such a numerical calculation of the radiative 
tail without approximations is presented in section In the present section, we give an 
analytical result for the soft-photon (i.e. small AEg) integral of Eq. ( |A43| ) , without any 



further approximations (Remark that in only an approximate evaluation of Eq. (|A43|) 
has been given). 

As is immediately seen by power counting, the integral in Eq. ( [A43| ) has a logarithmic IR 
divergence, corresponding with the emission of photons with zero energy. To demonstrate 
the cancellation with the IR divergence of the vertex diagram as stated above, one has to 
regularize the integral of Eq. ([A43|) . In this work this is performed by also using dimensional 



regularization. The soft photon integral is then evaluated in D — 1 dimensions {D ^ A 
corresponds to the physical limit). One now sees that it is extremely advantageous to have 
a spherical integration volume, in order to evaluate the integral for dimensions D ^ 4. 
Before continuing the integral of Eq. (|A43|) into D — 1 dimensions, the integration limits 



for / have to be made dimensionless, which leads in the dimensional regularization scheme 
to the introduction of the same scale /i in Eq. ( |A50|) as was introduced when changing 
the dimension of the virtual photon loop integral of Eq. (|A5|). This leads then in D — 1 
dimensions, to the bremsstrahlung integral : 
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(A50) 



The integral in Eq. ( |A50| ) is worked out by introducing polar coordinates in D — 1 
dimensions. To define the polar angle in the interference term of Eq. ( |A50|) , a Feynman 
parametrization is performed. This leads for I to the expression : 



•'\</^Es/^l jD-2 
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, (A51) 



1-P'^.l^ 

where / is the unit-vector along the soft photon direction, (3e = \Pe\, Pe = IPel the 

incoming and outgoing electron velocities (in the system S) respectively and where (3y = \f3y\ 
with 



I = — 



8' = ^ 



(A52) 



The integrals over 1 and the azimuthal angular integral (over D — 2 dimensions) can be 
performed immediately which yields : 
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The IR divergent term and the finite term are obtained by developing the polar angular 
integral in Eq. ([A53|) as 
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Performing the integrations in Eq. ( |A54| ) (the second integral in Eq. ( [A54| ) is simplified by 
making the substitution x— *>« = /?/(! — [3x) ) yields 



dx 



1 



{l-(3xY 



l-/?2 
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Consequently, the IR divergent term and the finite term of the integral / are obtained by 
using Eq. ( |A55|) in Eq. ( |A53| ) and by developing all other factors also to order ejR : 
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where the remaining Feynman parameter integrals ly^^ and ly'^' are given by 
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and where f3y is given by Eq. ( [A52| ). The integral Jy^^ in Eq. ([A57|) can be performed easily 
and yields 
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(A58) 



with V as defined in Eq. ( |A12| ). To obtain an anal5d;ical formula for the integral I^^ is much 
harder but was performed in Ref. E^, which we checked ^ and which yields the result : 
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(A59) 



where a, 0i and 02 are given by : 



cosh a 



l3e/3'^ sin 6e 



a > 0) 



Note that the relevant formula quoted in Ref. |47] contains some typing errors 
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cosh 01 = j3e cosh a , sinh 0i 

cosh 02 = cosh a , sinh 02 

The function $ in Eq. ( [A59| ) is given by 



-McOS^e+/3e' 

/3e/3g sin 6e 

PJ'^ cos fi^e - 



f3e(3'^ sin 6'e 



(A60) 
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which agrees with the Spence function (Eq. (|A13|) ) when x < 1. Compared with previous 
calculations in the literature, it was shown in Ref. that this integral ly"^^ was approxi- 
mated in Ref. and that the calculation of this integral in Ref. |^ contains a factor two 
error. We also checked the analytical formula of Eq. (|A59| ) by performing the integral of 
Eq. ( [A57| ) numerically. 

In the ultrarelativistic limit (/3e,/3l — > 1), the integral ly^^ of Eq. ( [A59| ) reduces to 
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Putting all pieces together, the result for the bremsstrahlung cross section accompanying 
elastic electron scattering is obtained as 
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where the finite part of the real radiative corrections is given by 



(A63) 
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where the expression of Eq. ( |A65|) in the » m? hmit has been rewritten in Eq. (|A66|) 



to allow comparison with other expressions found in the literature. 

Finally to evaluate we have to express the quantities in the system S in terms of lah 
quantities. The relations given in Eqs. (|A44|, [A45| )) yield for elastic scattering : = E'^, 



E'^ = Ee, and cos6'e = cos6'e. From the formula for 6r (e.g. Eq. (|A66|) , one then sees that 
one formally obtains exactly the same expression in terms of the lab quantities Ee,E'^,6e- 
The quantity AEg is calculated from the cut in the electron spectrum, using the expression 
of Eq. (|M9|) . 

A comparison of expressions Eqs. ( |A63|) , (|A64|) with the literature, shows that the same 
result is obtained as in Ref. A comparison with the expression used by Mo and Tsai 



22 1 will be given in the next section when we add the vertex correction and soft photon 



emission contribution, because only their sum is IR finite (and thus independent of the IR 
regularization procedure used). 



6. Elastic lepton-nucleon scattering 



In this section, we bring together the first order radiative corrections at the lepton 
side (lepton vertex and soft bremsstrahlung from the lepton) and the photon polarization 
correction to correct the elastic lepton-nucleon scattering cross section. As was shown in 
the previous sections, these corrections can be calculated model-independently. In the next 
section, we discuss the additional radiative corrections to the lepton-proton cross section, 
which originate from the proton side (proton vertex correction, soft bremsstrahlung from 
proton and two-photon exchange corrections). To calculate these corrections at the proton 
side, a model for the off-shell (or half off-shell) 'jNN vertex is needed however, and which is 
therefore to some extent model-dependent. For this latter part, we will refer to the recent 
work of Ref. H. 

The elastic lepton scattering cross section, corrected to first order in aem for the lepton 
vertex contribution and for the photon polarization contribution, is given by 
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where the finite part Syertex of the lepton vertex correction is found from Eq. ( [A17|) to be 
given by 
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(A68) 



In writing down Eq. ( [A67| ) to first order in Oem; the contribution of the anomalous mag- 



netic moment term G (Q^) in the vertex correction Eq. i \KE[ ] has been dropped. This contri 
bution vanishes in the ultrarelativistic hmit { » ) as can be seen from Eq. ( All ) 



The first term in the last line of Eq. ( [A68|) corresponds with the vertex correction term 
quoted by Mo and Tsai (Eq. ( II.5) of Ref. H). 

The finite part of the photon polarization correction, 5^ac = 2n((5^), follows from 
Eq. ( [A40| ) as 
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which agrees with the expression quoted by Mo and Tsai (Eq. ( II. 4) of Ref. |2^). To 
evaluate the vacuum polarization due to fj^~^fi~ and t~^t~ pairs at intermediate Q', one has 
to use Eq. (|A69D instead of the limit of Eq. (|A70|) . Q 

When adding the real (Eq. ( |A63| )) and virtual (Eq. ( |A67| )) radiative corrections at the 
lepton side, one verifies that the IR divergent parts exactly cancel. The remaining finite 
contribution is given to first order in aem by 
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where 6yac, ^vertex and are given by Eqs. ( [A69| ), ( |A68D , and ( |A64| )-( [A^ respectively. 
Bringing the three contributions together, leads to the expression (in the Q' » m? limit) 
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where AEg, which is the maximum soft-photon energy in the cm. system of (recoiling 
proton + soft-photon), is determined as in Eq. ( |A49| ) , when applying this formula to the 



^Note that an incorrect expression is used in |50| for the vacuum polarization contribution due to 
pairs (Eq. (A5) in their paper). 
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scattered electron spectrum. 

We can compare Eq. ( A72| ) with the recent calculation of Maximon and Tjon where 
this calculation was also performed (using a finite photon mass to regularize the IR diver- 
gences) without doing any approximations. Comparing Eq. (|A72|) with their Z-independent 
term [Z being the hadron charge) - i.e. when not considering radiative corrections at the 
proton side or two-photon exchange contributions at this point - we find exactly the same 
result. As was noted in Ref. pi], the last two terms of Eq. ([A72|) were omitted by Mo and 



Tsai \22 



We can approximately take into account the higher order radiative corrections by ex- 
ponentiating the first order vertex and real radiative corrections. This is strictly true only 
for the IR divergent part of the vertex correction and soft photon emission contribution, 

mm 



and was demonstrated in Refs. 



(see e.g. Refs. |51|J5^ for pedagogical derivations). 
The application of this exponentiation procedure also to the finite part consists of an ap- 
proximation which can be checked by comparing the result with the first order formula 
of Eq. (|A71|) . For the photon polarization contribution, we iterate the first order vacuum 
polarization contribution of Eq. ( [A69| ) to all orders (resumming all vacuum bubbles of the 
type of Fig. |T9| (c)) by keeping the photon self-energy in the denominator as in Eq. ([A67D . 
Remark that a resummation of the first order vacuum polarization contribution does not 
lead to an exponentiated form. Assuming exponentiation for the finite parts of the vertex 
and soft photon emission contributions - as occurs for their IR divergent pieces - leads then 
to the radiative correction formula 



d(j 



dn: 



+ 



e/ VIRTUAL-f 



da 



dn: 



e/ REAL SOFT ^ 



da 

dnL 



J ^vertex ~h 



BORN 



(A73) 



7. Radiative corrections at the proton side and two-photon exchange contributions 



In the previous sections, we considered radiative corrections to elastic electron scattering 
originating solely from the electron side (vertex correction and bremsstrahlung) and from 
the vacuum polarization. These corrections, which are the dominant ones, can be calculated 
model independently and follow from QED. To calculate the first order radiative corrections 
originating from the proton side (proton vertex correction, bremsstrahlung from proton and 
direct and crossed two-photon exchange contributions), one needs a model for the internal 
structure of the nucleon because one requires knowledge of off-shell (or half off-shell) 'jNN 
vertices. This model dependence will become important if one aims at a precision of electron 
scattering experiments at the 1 % level. To quantify the magnitude of those effects, we refer 
to the recent work of Maximon and Tjon [^], where an initial study was performed of the 
size of internal structure effects. 

In Ref. ^ 



the proton current was taken to have the usual on-shell form and form 



factors were included in the calculation. The calculation of Ref. ||2J] goes beyond previous 
works [^,^], as the proton vertex correction and the bremsstrahlung from the proton where 
calculated without approximations within the given model for the proton current. In the 
calculation of the direct and crossed box diagrams (two-photon exchange contributions), a 
less drastic approximation was made in [|2^ as in [^] (where those box diagrams where only 



calculated in the soft-photon approximation). 
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The calculation of Ref. yields then the correction formula for elastic electron scat- 
tering : 
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where 6yac, Syertex and 6r are given as above (Eg. ( [A7^ )). The terms in Eq. ( |A74| ) proportional 
to Z (hadron charge) and contain the corrections from the proton side. The correction 
6i, proportional to Z, was calculated in Ref. 



as 
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where lS.Es and ?7 are given as in Eq. (|A49|) and where the variable x is defined by 
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The correction proportional to Z^ was split into two parts in Ref. The contribution 

independent of the nucleon form factors was calculated in Ref. as : 
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where p is defined as in Eq. ( A76 ), and where E'j^ {\p'n\) ^^e lah energy (momentum) 
of the recoiling nucleon. For the lengthier expression of which depends on the nucleon 



form factors, we refer to Ref. 24 
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APPENDIX B: TREATMENT OF SINGULARITIES 



In the numerical calculation of the amplitudes for the virtual photon radiative correc- 
tions to the ep — >■ ep7 reaction, we need to calculate two or three dimensional Feynman 
parameter integrals, as discussed in section [III Q . In the integration over the first variable, 
the numerator consists of polynomials and the denominators may have some structures of 
the form {a'x+(3'±ie')"', or (a'x^+/?'x+7'±2£:')" with n=l,2. Therefore, in the calculations, 
the following integrals appear : 



lim 



x™ dx 



(a'x + f3'± ie') 



or 



lim 

£'^0 + 



x"^ dx 



{a'x^ + (5'x + 7' ± ie') 



fBll 



When the denominator has no singularities in the integration range, it is, in principle, easy 
to calculate these integrals which have the form 



x"^ dx 
(a'x + 



or 



x™' dx 



„ (a'x2 + fi'x + 7')" ' 



(B2) 



Some recurrence relations for these integrals are known [ |5B| , |5^ , but for small values of a' as 
compared to or to 7', these relations are unstable numerically. In these cases, we have 
used either a Taylor expansion or the usual Gauss- Legendre integration method to get very 
accurate results. 

In the following part of this appendix, we give the relations used when the denominators 
in Eq. ( |B1[ ) have singularities in the integration range except in a or h. The details are given 
elsewhere ||28|. The principle of the method is based on the following relation : 
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Each integral can be separated in a real part and an imaginary part and we can use for them 
the analytical expressions given in . 

Let us start with the case where the denominator is a polynomial of degree 1 in the 
integration variable. In that case, there is only one singularity for xo = —(3' /a' and the sign 
of the imaginary part will depend on the sign of a' . For n = 1 and a' > 0, we have 



lim 
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When a' < 0, we have only to replace ±ze by ~^ie in the right hand side of the Eq. (]B4| 
We now define the following quantities 
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n > 2 



to obtain the relations 
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For n = 2 and a' > 0, we have 
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When a' < 0, we have only to replace ±z£: by ^ie in the right hand side of the Eq. ( |B11D . 
We next define the following quantities : 
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In terms of these quantities, the integrals of Eq. ( |B1| ) with n = 2 are given by 
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We can notice that the real part of these integrals for n = 1 as well as for n = 2 can be 
derived from the binomial expansion (xq + X)™. In the case of n = 1, the imaginary part 
is proportional to 7r/(xo) where f{x) is the numerator of the integrand. For n = 2, it is 
straightforward to show p8[ that the imaginary part is proportional to 7r/'(xo). 

When the form of the denominator is (a'x^ + + 7' ± ie')^, i.e. a polynomial of degree 
2 in the integration variable, it is always possible to come back to the preceding cases. When 
a' > 0, we have 

x"^ dx 1 x"^ dx , 
hill / -, ^ = — - lim / — — ^ , (B19) 

e'^o+ {a'x^ + 13' X + 7' ± le'Y a'" e^o+ J ^ {x^ + I3x + ±ieY ^ ' 

with the following definitions : 

3' 7' e' 

(3=-,. 7 = ^, e = -^. (B20) 
a' a' a' 

The integrand in Eq. ( |B19D has some singularities when 5 = — 47 is positive. When 
a' < 0, we have only to replace ±i6 by ^ie on the right hand side of Eq. (|B19| ). 
It can be shown |]2S| that 

X™ rfx /■'' x"" dx 
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These integrals can be easily calculated using the decomposition of the fraction into elemen- 
tary fractions. For n = 1, we obtain 

dx 
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and for n = 2 
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TABLES 



TABLE I. Radiative corrections to elastic electron-proton scattering for MAMI and JLab kine- 
matics. First column : in GcV, second column : Of, in deg, third column : in GeV^. See 
text for details on the different contributions. The real radiative corrections are calculated with 
[E'^'' — E'^) = 0.01 E^. The total radiative correction (to first order) is indicated by 5tot, and the 
exponentiated (EXP) result (except for the vacuum polarization contribution, see text) is shown 
in the last column. 
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FIGURES 




(BORNi) (BORNf) 



FIG. 1. Tree level contributions to the ep — > epj reaction: Bethe-Heitler diagrams (a) and (b); 
nucleon Born diagrams (c) and (d). 
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FIG. 2. First order virtual photon radiative corrections to the ep ejrf reaction. 
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FIG. 4. Test of convergence : The relative effect of the virtual radiative corrections to the 
BH + Born cross section is shown as a function of the angle between the two photons {q, q') for 
a typical MAMI kinematics. The curves correspond to tests performed with different densities 
of integration zones and points near the edge of the integration domain. They show the good 
numerical convergence obtained. 
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FIG. 5. Analogous test of convergence as in Fig. ^ but for JLab kinematics. 
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FIG. 6. Test of gauge invariance for MAMI kinematics. The dotted curve shows (for ihustrative 
purpose only) the result when only the diagrams (V2i) and (V2f ) of Fig. ^ are included. The dashed 
curve is the result of all analytically calculated virtual radiative corrections. The full curves show 
the result when also the numerical contributions (Feynman parameter integrals) are included. 
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Evolution in mass 
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FIG. 7. Mass dependence of the virtual radiative corrections for MAMI kinematics. 



63 




FIG. 8. Detail of the cross section for photon emission from an electron (Bethe-Heitler cross sec- 
tion), when the photon is emitted around the incoming electron direction. The electron kinematics 
correspond with : Ee = 855.0 MeV, E'^ = 621.4 MeV, Oe = 52.18°. 
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elastic e-p scattering: £^=855.00 MeV, i3e=52.180° 




FIG. 9. Radiative tail of elastic electron-proton scattering. Upper plot : fully numerical cal- 
culation (black points) compared with the 1/ASg dependence of the soft photon result (straight 
line). Lower plot : deviation between the full calculation, when only radiation from the electron 
is included (open diamonds) and when both radiation from electron and proton are taken into 
account (black points), with the soft photon result (straight lines). See text for details. 
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elastic e-p scattering: Ee=855.00 MeV, de=52.180° 
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FIG. 10. Radiative tail of elastic electron-proton scattering at Eg = 855.00 MeV and 9e = 
52.18°. A comparison is shown between fully numerical calculation (indicated by the points) and 
the simulation (curves, see text). 
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FIG. 12. Determination of the elastic cross section for the kinematics Ee = 705.11 MeV and 9e 
42.6°. 
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FIG. 13. Differential ep ejry cross section for MAMI kinematics at g =33 MeV/c. 
Dashed-dotted curve : BH + Born contribution, dashed curve : BH + Born + virtual radia- 
tive correction, full curve : BH + Born + total radiative correction. The real radiative correction 
is shown here for a maximal soft-photon energy of A£?s = 10 MeV. 
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FIG. 14. Differential ep epj cross section for MAMI kinematics at q = 111.5 MeV/c. Curve 
conventions as in Fig.^. The real radiative correction is shown here for a maximal soft-photon 
energy of AEg = 10 MeV. 
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FIG. 15. Differential ep ep^ cross section for JLab kinematics at q' = 120 MeV/c. The BH 
+ Born result is compared with the result including virtual radiative corrections. 
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q = 600 MeV/c, q' = 120 MeV/c, 8 = 0.61 
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FIG. 16. Double polarization asymmetry for VCS with proton polarized along the z-axis (upper 
panel) or polarized along the x-axis (lower panel) for MAMI kinematics. Dashed-dotted curve : 
BH -I- Born, full curve : BH -|- Born -|- radiative corrections (both curves nearly coincide!). 
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FIG. 17. Differential ep ep-f cross section in lab : DVCS kinematics at JLab. Dotted curve 
: BH, dashed curve : DVCS, dashed-dotted curve : BH + DVCS, thin full curve : BH + DVCS 
+ virtual radiative corrections. The thick full curve represents the BH + DVCS + virtual and 
real radiative corrections, where the real radiative corrections are calculated with AEg = 0.1 GeV, 
which corresponds with a cut in the missing mass spectrum (Eq. (|59[)) of M^^^ - ~ 0.21 GeV^. 
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FIG. 18. Electron single spin asymmetry : DVCS kinematics at JLab. Dashed-dotted 
BH + DVCS, full curve : BH + DVCS + radiative corrections. 
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FIG. 19. First order virtual and real radiative correction processes : (a) vertex diagram, (b) 
lepton self energy diagram, (c) photon polarization diagram, (d) and (e) soft-photon emission 
contributions to elastic lepton-nucleon scattering. 
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